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We calculate the tunnelling current through a Fabry-Perot interferometer in the fractional quan- 
tum Hall regime. Within linear response theory (weak tunnelling but arbitrary source-drain voltage) 
we find a general expression for the current due to tunnelling of quasiparticles in terms of Carlson's 
R function. Our result is valid for fractional quantum Hall states with an edge theory consisting of 
a charged channel and any number of neutral channels, with possibly different edge velocities and 
different chiralities. We analyse the case with a single neutral channel in detail, which applies for 
instance to the edge of the Moore-Read state. In addition we consider an asymmetric interferom- 
eter with different edge lengths between the point contacts on opposite edges, and we study the 
behaviour of the current as a function of varying edge length. Recent experiments attempted to 
measure the Aharanov-Bohm effect by changing the area inside the interferometer using a plunger 
gate. Theoretical analyses of these experiments have so far not taken into account the accompanying 
change in the edge lengths. We show that the tunnelling current exhibits multiple osculations as a 
function of this edge length, with frequencies proportional to the injected edge current and inversely 
proportional to the edge velocities. In particular the edge velocities can be measured by looking at 
the Fourier spectrum of the edge current. We provide a numerical scheme to calculate and plot the 
R function, and include sample plots for a variety of edge states with parameter values which are 
experimentally relevant. 



I. INTRODUCTION 



The physics of the fractional quantum Hall effect is not 
captured through the conventional picture of symmetry 
breaking and local order parameters. The effect arises 
in low-dimensional electronic systems and is a prime ex- 
ample of a topological phase of matter^"^. A range of 
fractional quantum Hall phases has been discovered''"^, 
each characterized by its filling fraction v, which deter- 
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mines the Hall conductivity through an ~ ^^fi- Each 
of these phases possess a different type of order known 
as topological order*^, referring to the presence of long- 
range entanglement of the ground state'^. Some manifes- 
tations of topological order in the FQHE are the simulta- 
neous formation of an energy gap in the bulk and gapless 
states at the edge of the system^"-'^, a ground-state de- 
generacy determined by the topology of the space-time 
manifold^^'^'^ and exotic properties of the low-lying exci- 
tations of the system^"*''^. 

In particular, topological order predicts quasiparticle 
excitations known as anyons'^'^^, which possess frac- 
tional charge^'*'^*''^ and obey a generalized form of 
exchange statistics^^'^°'^^. These statistics generalize 
bosonic and fermionic statistics in the sense that inter- 
change of anyons multiplies the wavefunction by a phase 
factor (Abelian anyons) or induces a rotation in an in- 



ternal, non-local space of degenerate states (non- Abelian 
anyons). Non- Abelian anyons have been put forward as 
candidates for the realization of a topological quantum 
computer^'^^"^"'. 

A candidate for a non- Abelian FQH phase is the exper- 
imentally observed v = \ state^'' . The Moore- Read Pfaf- 
fian state^^'^^ and its particle-hole conjugate the Anti- 
Pfaffian^^''^^ could possible describe the corresponding 
topological order. These states predict the same Hall 
conductivity, but differ in other properties such as the 
type of anyons present and the effective edge theory of 
the system. It is both a theoretical and experimental 
challenge to design experiments sensitive to the topolog- 
ical order of the system, thereby identifying the nature 
of the I state and other potential non- Abelian phases. 

Experiments have successfully measured the fractional 
charge of tunnelling quasiparticles'®'"'^^'^^ for a variety 
of quantum Hall phases. More recent experiments aim 
to fully determine the topological order through use of 
various interferometric devices^^""^^ . These experiments 
make use of the transport properties through these de- 
vices, which is determined by the edge where the electric 
current is located'^'''^. The electric current is chiral and 
flows along the edge in a single direction. Backscattering 
can occur through quantum point contacts, where oppo- 
site edges are forced together. This induces tunnelling of 
anyons between the edges^®''^^. The resulting tunnelling 






Figure 1. Figure of an interferometer. Constrictions bring 
the edges together and causes tunnelling of charge from one 
edge to the other. Tunnelling occurs only between the inner 
edge states. The outer edges carry excitations of lower density 
Hall liquids, such as underlying fully filled Landau levels. We 
assume these are fully transmitted. 



or backscattering current through the constrictions, Ib, 
depends on the type of anyon tunncUing. 

In this work we analyse the tunnelhng current 
through a Fabry-Perot interferometer in linear response 
theory "^^"^^ . An interferometer consists of multiple points 
contacts, see Figure I. Anyons tunnel along different tra- 
jectories, which gives rise to interference effects. In a 
simple picture we assign ti and ^2 as the complex ampli- 
tude of a quasiparticlc tunnelling along the correspond- 
ing point contact. The tunnelling current follows from 
the absolute value 



Ib - \ti + 12|' - \tif + \t2f + 2Re[iit2] 



(1) 



In linear response theory the form |iip and |i2p is radi- 
cally different for the case of a tunnelhng anyon as com- 
pared to what would be expected for electrons. It is a 
non-linear function of the applied voltage, the tempera- 
ture of the system, and the fractional charge and scal- 
ing dimension of the tunnelling anyon'^*''*^. The term 
2Re[iit2] is the interference current. Interference arises 
due to a variety of causes, such as the Aharonov-Bohm 
effect, the relative phases of the tunnelling coupling con- 
stants and the dynamical interference due to the finite 
velocity of the anyons traversing the interferometer. Per- 
haps the most interesting contribution to the interference 
current is due to the statistics of the anyons. Anyons lo- 
calized in the bulk and inside the interferometer braid 
with anyons tunnelling between the edges. This braid- 
ing of anyons effectively reads out the topological state 
of the bulk anyons, and this signature manifests itself in 
the interference current '*^"^^. Further effects arise that 
go beyond braiding properties which are due to coupling 
of bulk quasiparticles and edge degrees of freedom'"''^"^^. 

We are primarily interested in the dependence of the 
interference current on the dynamical properties of the 
edge, such as the velocity of the edge channels and the 
applied voltage. Earlier work focused on edge states with 
a single characteristic velocity ^° or edge states of specific 
quantum Hall candidates ' to obtain an expression for 
the interference current. We present here the more gen- 
eral case of an asymmetric interferometer, a generic num- 



ber of edge channels with possibly different edge veloci- 
ties and opposite chiralities, at both zero and finite tem- 
peratures. 

Our result is an analytic expression for the interference 
current in terms of a generalized hypergeometric func- 
tion known as Carlson's R function^"^. This scaling func- 
tion is closely related to the Lauricella hypergeometric 
function^^'^^ . This Lauricella function is a multivariable 
generalization of the Gauss hypergeometric function^°, 
which is a function which enters the expression for the 
interference current for edge states described by a sin- 
gle velocity *°. Our expression generalizes this result for 
edge states consisting of an arbitrary number of decou- 
pled channels described in the conformal limit. Each of 
these channels has its own corresponding velocity. We 
also find an expression for the interference current at zero 
temperature in terms of the confluent Lauricella hyperge- 
ometric function^^ , which is a multivariable generaliza- 
tion of the Bessel function of the first kind. Finally we 
obtain an expression for the two-point correlator of an 
anyon situated at the edge in the (a;,x)-representation. 

As a function of the voltage between the two edges the 
interference current behaves as a sum of decaying oscil- 
lations. The frequencies of these oscillations are deter- 
mined by the edge lengths, edge velocities and the quasi- 
particlc charge. For an antisymmetric interferometer this 
results in four frequencies appearing in the Fourier spec- 
trum of the interference current as a function of the volt- 
age. These four frequencies correspond to the possible 
combinations of one edge length and one edge velocity. 
Alternatively, we can fix the voltage and vary the length 
of one edge. This again results in oscillating behaviour 
with frequencies determined by the voltage, edge veloci- 
ties and the quasiparticle charge. 

This behaviour of the interference current as a func- 
tion of varying the edge length is relevant to experiments 
which measure the Aharonov-Bohm oscillations through 
application of a plunger gatc'^^'^^"'^^. The plunger gate 
effectively deforms the area inside the interferometer 
through use of the Coulomb interaction. This deforma- 
tion changes the Aharonov-Bohm phase of the tunnelling 
quasiparticles, which results in an oscillating interference 
current as a function of the side-gate voltage. The fre- 
quency of these oscillations, which we denote by 4>ab, can 
be used to measure the charge of the tunnelling quasipar- 
ticle and the effect of quasiparticle braiding*^''*'^"'''^. 

However, the change in area of the interferometer can 
also result in a change in the edge length, depending 
on the specific geometry of the interferometer. We show 
that for certain assumptions, such as the geometry of the 
device, this change in edge length results in additional 
oscillations in the interference current. When the change 
in edge length is large enough and linear with the side- 
gate voltage, then the interference current shows multiple 
oscillations characterized by the frequencies ip^B , 'Pab + 
^^ and (JJAB + ^rii- These shifted frequencies can be 
used to measure the edge velocity. 

The paper is structured as follows. We start in Sec- 



tion II with a discussion of the edge theory of a fractional 
quantum Hall phase. We specify the structure of the edge 
theory and quasiparticle operators, which is based on the 
decomposition in terms of a charged and neutral channel. 

In Section III we discuss the model Hamiltonian of the 
Fabry-Perot interferometer in terms of the quasiparticle 
operators and the corresponding linear response. This 
leads to an expression of the tunnelling current in terms 
of four-point correlators of the quasiparticle operators, as 
shown in Section IV. Specifically, the tunnelling current 
is given by evaluating the Fourier transform of these four- 
point correlators at the value of the Josephson frequency, 
see expression (41). 

The four-point correlators depend on the precise na- 
ture of the edge state and they do not have a universal 
form. But as we show in Section V the correlator has a 
leading dependence which does have a universal expres- 
sion, which is a result of the conformal symmetry in the 
large system-size limit. 

This leads to our main result in Section VI, which is 
the Fourier transform of the leading order expression of 
the four-point correlators at finite temperature, Eq. (68). 
This expression is given in terms of Carlson's R function 
which acts as a modulating function. Since this func- 
tion is somewhat obscure we summarize its properties in 
Appendix B3 and describe our method of computing the 
function, which is through its relation to the Lauricella 
function. 

In Section VII the main result is further explored for 
special cases, such as the zero temperature case. In Sec- 
tion VIII we plot the interference current and the R func- 
tion for a range of experimentally relevant parameters 
and analyse the result for a number of trial states for the 
ly = 5/2, v = 7/3 and i^ = 12/5 plateaus. In general 
the R function has a decaying oscillating behaviour. We 
show how the frequencies of these oscillations relate to 
the physical parameters of the system. In Section IX we 
discuss the relevance of our results to experiments involv- 
ing the Aharonov-Bohm phase in the interferometer. 



II. EDGE THEORY 

A quantum Hall fluid is an example of a topological 
system^""^'^. The fluid has a mobility gap in the bulk of 
the system. Simultaneously, gapless states develop at the 
edge of the system where the confining potential crosses 
the Fermi level^°'^^. These gapless edge states are chiral 
and responsible for the transport properties of the fluid. 

The effective edge theory of the fractional quantum 
Hall effect can be seen as as a consequence of anomaly 
cancellation^^'^^'^^'^^"^^. The effective bulk theory of a 
quantum Hall fluid is a Chern-Simons theory; a topo- 
logical gauge theory which describes the bulk of the sys- 
tem and develops an anomaly on the boundary where 
gauge invariance is broken. A dynamical edge theory 
forms, with the same anomaly, but opposite in sign. The 
combined bulk plus edge system is gauge invariant and 



anomaly free. 

In the long-wavelength approximation the resulting 
edge theory is a chiral conformal held theory. The elec- 
tron and quasiparticles of the theory are represented by 
local operators in this conformal fleld theory. The set of 
all local operators forms the chiral algebra®^. By speci- 
fying the chiral algebra we zoom in on a candidate frac- 
tional quantum Hall state at some flUing fraction i'. To 
be a suitable candidate for a quantum Hall state, the chi- 
ral algebra needs to fulfll a number of conditions. These 
conditions include for instance the existence of an elec- 
tron operator and the presence of a U{1) symmetry. We 
assume such conditions are always satisfled in our discus- 
sion. 

The U{1) symmetry arises due to presence of the elec- 
tric current. In the case of a Laughlin state the U{1) 
symmetry is the full gauge symmetry of the bulk Chern- 
Simons theory. The corresponding edge theory is a chi- 
ral u(l) current algebra, also known as the chiral bo- 
son or chiral Luttinger liquid^^'^*"^°. More complicated 
Abelian states involve the presence of multiple chiral 
bosons^'^^'^^'^^. For non- Abelian quantum Hall states 
the f/(l) symmetry is also present, but only as a subgroup 
of a larger, more complicated gauge group^'^''*^'^^'^'''^"'. 
Following Ref. 65 and 72 we limit ourselves to those 
states described by a representation of an algebra which 
is formed by a direct product 



A^Wn(^ U{1) 



(2) 



Here yV„ is the symmetry of the chiral algebra responsible 
for the non- Abelian nature of the system. Quasiparticle 
operators obey the same decomposition. We refer to the 
different terms in the product as the neutral and charged 
channel of the edge theory. Throughout the main text 
we mostly deal with a single charged and a single neutral 
channel, although we comment on the more general case 
of edge states with multiple modes. 

Frequently, we deal with quantum Hall states which 
develop on top of one or multiple completely flUed Lan- 
dau levels, as is the case with candidate states for the 
fllling fraction i/ = 5/2. These filled Landau levels form 
edge states as well, and for simplicity we assume these 
states completely decouple from the quantum Halls state 
of interest. In the presence of a point contact these filled 
edge states are assumed to fully transmit, meaning charge 
transfers only between the inner-most edges, see Figure I. 



A. Charged channel — the chiral boson 

The action of the charged channel is that of the chiral 
boson'^^'*^'^^"™. We consider a single edge with a right- 
moving chiral boson, held at a voltage bias U in the gauge 



ttx — 0. The action is given by 

S'c = — / dtdx [dttpd^tp - Vc{dxtpf 



I — /• 

eU / dtdx \dx'-p\ ■ 

2tt J 



(3) 



The field is compactified by the identification 
Lp — if -\- 2ttv and f c is the velocity of the channel. 
The field ip represents the charge density along the edge 
through the relation 



p{x) 



2lT 



dx^ 



(4) 



Quantization^^ results in the (non-local) equal-time com- 
mutation relations 



[ip{x), ^{y)] = -iT:sgn{x - y) 
[dxip{x), ip{y)] = -i27r(5(a; - y) 



(5) 



with sgn(a;) = +1,0,-1 for x > 0, a; = and x < 
respectively. Hamilton's equations of motion are given 

by 



{dt - Vcdx)(p = -^/veU . 



(6) 



Differentiating with respect to x shows the charge den- 
sity p is a conserved current as long as the edge is an 
equipotential, [dt — Vcdx)p = 0. The corresponding con- 
served charge is (up to normalization) identified as the 
total charge operator 



p{x)dx 



27r 



dx'p{x)dx 



(7) 



The Hamiltonian Ki^ q,c for a right moving edge held sub- 
ject to the potential U which follows from the action (3) 
is 



K 



LM,c 



= H^ 



L,0,c 



-eUQ 



(8) 



The second term, eUQ, is the coupling to the electro- 
static potential. The first term corresponds to the Hamil- 
tonian of the system in the absence of an external poten- 
tial, 



H 



L,0,c 



An 



dx{dxvf 



(9) 



The full Hamiltonian (8) is a generalization of the usual 
grand canonical Hamiltonian of the form Kq = Hq — /xiV, 
with N the number operator. Instead of a number oper- 
ator we use the charge operator. 



B. Neutral channel and quasiparticles 

We do not explicitly specify the nature of the neutral 
channel, but only assume the decomposition (2). What 



matters is that the full chiral algebra fixes the quasi- 
particle content of the theory and it comes equipped 
with consistent rules for fusion and braiding of these 
quasiparticles^'^'^. Each quasiparticle is characterized by 
its conformal dimension and its fusion and braiding rules 
with respect to the remaining quasiparticles. This spec- 
ifies its quantum dimension as well. 

A general quasiparticle operator factorizes as 



■0 {x, t) (X a{x, t) ® e 



^(p{x^t) 



(10) 



The exponentiated operator e'^''' is normal ordered 
and corresponds to the charged channel, while a repre- 
sents the neutral channel. The normalization factor de- 
pends on the regulators of the theory^" . The operator is 
characterized by its conformal dimension, /i^t = h„ + he- 
The conformal dimension of the charged channel follows 
from the charge and the filling fraction, hc= ^- The 
commutation relations (5) show that the operator obeys 



[Q,ijHx,t)]^Qi;^x,t) 



(11) 



and so the corresponding quasiparticle carries an electric 
charge Qe. 

For each quasiparticle a conjugate particle exists with 
opposite charge and the same conformal dimension^"'. We 
set 



^{x,t) — (T{x,t) (g) e 



^^{x,t) 



(12) 



The operator a is chosen such that the fusion product of 
a and a contains the identity channel, 

(jxa^l + ... . (13) 

For non-Abelian quasiparticles we have, in general, mul- 
tiple fusion channels. We assume that, for a given neutral 
mode, for each operator a there is a unique conjugate op- 
erator a in the theory which obeys the fusion rule (13). 
This assumption is in fact a condition on the chiral alge- 
bra. 

Finally, we also mention that the neutral channel 
traverses at some characteristic velocity w„ and it is 
equipped with some neutral Hamiltonian, iJ„, similar to 
the charged channel. However, the neutral channel does 
not couple to the electromagnetic field, and therefore no 
analogous coupling of a zero mode to the external elec- 
trostatic potential appears. Furthermore we assume the 
general situation in which u„ 7^ Uc. 



III. MODEL OF A FABRY-PEROT 
INTERFEROMETER 

A. Tunnelling Hamiltonian 

In this section we treat the basic idea behind the tun- 
nelling formalism in a system of point contacts'*"'''^'''^'^'^. 
We consider a quantum Hall bar of infinite length at a 




Figure 2. Figure of an interferometer. Tunnelling of quasipar- 
ticles occurs at the point contacts, e.g. from Xi to yi through 
the operator Vi. The dotted arrows represent the direction of 
the edge currents, with a right moving current on the lower 
In the text we set a = \yj — yi\ and b — \xj — a;^ | . 



Schroedinger picture for the tunnelling Hamiltonian 

Ht = T + T^ (18) 

where the T operator is defined as 

N 

T = ^r,V(x„yO ■ (19) 

i=l 

Here the sum runs over the N point contacts and xi and 
yi denote the coordinate of the z'th point contact on the 
lower and upper edge. 

B. Tunnelling Current 



uniform filling fi'action v. The two edges, denoted as 
E/j/i, are disconnected and multiple constrictions are de- 
scribed by hopping terms allowing for the tunnelling of 
quasiparticles from one edge to the other. Here the sub- 
script L and R denote the left (upper) and right moving 
(lower) edge of the system. For each edge we have an 
electric charge operator 



Qr,i 



dxpR/L{x) . 



(14) 



We apply a voltage bias between the two edges, which 
is incorporated by fixing the electrostatic potentials Uji 
and Ul at the lower and upper edge respectively. The 
full Hamiltonian K is given by 

K = Ko+Ht . (15) 

K„ = Ha - cUlQl - eUiiQji (16) 

Here H^ is the tunnelling Hamiltonian which is treated 
perturbatively with respect to Kq. The grand-canonical 
Hamiltonian Kq consists of the terms coupling to the 
DC voltages through the charge operators and Hq. The 
Hamiltonian Hq decomposes into the Hamiltonians for 
the decoupled left and right moving edges, Hl and Hjj. 
In addition Hj^/fj describes both the charged and neutral 
channels He and iJ„ of each edge. 

The tunnelling Hamiltonian Ht couples the edges 
through tunnelling of (quasi)holes and (quasi)electrons. 
For this we first introduce the tunnelling operators V. 
We set X and y as the coordinates of the lower and upper 
edge respectively. A generic operator which tunnels a 
quasiparticle with charge e* = Qe (e > 0) from the lower 
to the upper edge is then 



y 



V{x,y) = V^(y)V'(a;) 



(17) 



The operators ■0 and V'^ are related as explained in Sec- 
tion II B. Similarly V\x, y) — ip^{x)'ip{y) tunnels a quasi- 
particle from the upper to the lower edge. 

We now consider a system of N well-separated point 
contacts. Each point contact is approximated by a 
single tunnelling operator V{xi,yi) and a correspond- 
ing tunnelling coupling constant T^. We have in the 



The quantity of interest is the current running through 
the point contacts from one edge to the other, the so- 
called backscattering or tunnelling current (/s). It is 
defined as the rate of change of the difference in electric 
charge of the edges, |J^ {Qr — Ql)- Using the equations 
of motion for operators in the Heisenberg picture we have 



If 



4[Qh 



QL,T + Tn 



(20) 



Here we used that the charge operators commute with the 
free Hamiltonian Hq as the charge is conserved separately 
on each edge in the unperturbed system. The commu- 
tation relations (11) imply [Qr,T] = -QT = - [Ql,T], 
and so we obtain 



iB = iQe{T-T^) 



C. Linear Response 



(21) 



Initially, at some reference time to, the perturbation 
Ht is absent and the two edges are decoupled. At this 
initial time to both edges are in thermal and (separate) 
chemical equilibrium with respect to the Hamiltonian Kq. 
The density matrix is given by 



Wo = w{to) 



^ e-^^° 



/Z 



(22) 



Note that the external DC voltage is not treated pertur- 
batively, but directly incorporated into the initial density 
matrix. 

The perturbation Ht is adiabatically switched on at 
t > to, slowly driving the system out of equilibrium. The 
time evolution follows from the usual time evolution op- 
erator UK{t,to) which solves the Schroedinger equation 
with respect to K, 



tdtUK{t,to) = KUK{t,to) . 



(23) 



In a perturbative approach^^ we factorize the time evo- 
lution operator as UK{t,to) = e-'^°^*-*o^UHT{t'){t,to)- 
Through (23) it follows that UnTif) satisfies 



ldtUH^(t'){tM - HT{t')UH^if){tM) 



Hrit) = e*-^«'i?Te 



-iKnt 



(24) 



The operator UHj.{t') is expanded as a power series in the 
tunnehing couphng constants. At first order in T{x, y) we 
have 



UHHt')it^to) = 1 - W dt'HTit') + ... (25) 

Jto 

The expectation value of an operator O is 
{0{t)) =Ti[woOK{t)] where Oxit) is the Heisen- 
berg representation of the operator 

OKit) - UJ,it,to)OK{tQ)UK{t,to) . (26) 

At the initial time to the perturbation Ht is absent, and 
so Ok {to) —OKo{to)- This identity together with the 
factorization of Uk and expression (2G) results in 

OK{t) = U%^^,,^(tM)OKo{t)UH^it'){tM) (27) 

where we have defined 

OK,,{t)=e'^°'Ose-''^''' (28) 

and Os is the Schroedinger picture of the operator. 
When UHT{t') is expanded and we keep only the lowest 
order term we obtain for the expectation value 



{0{t)) = {OKo{t)), 



dt'{[OKo{t),HT{t')]), 



(29) 



Here (• • •)p = Tr [wq • ■ ■ ] is the ensemble average with 
respect to the unperturbed thermal state of the edges, 
Eq. (22), and we have set to -> — cx). We emphasize that 
this thermal state still includes the nonperturbative effect 
of the DC voltage. Expression (29) is the Kubo formula 
for the operator O with respect to the perturbation Ht- 



where VHo{x,y,t) = e'^°*Vsix,y,t)e-'"''*. This is sim- 
plified further by using that when [A, B] = aB then 
g-»/3A^g»/3A ^ Be'"!^. This gives 



VKo{x,y,t) = e''^^'VH,{x,y,t). 



(31) 



Here we have defined ojq = Qe{Uii ^ Ul)/^, which is the 
Josephson frequency for a particle with charge Qe. The 
value of the charge Q depends on the specific edge and 
quasiparticle under consideration. Typical experiments 
are carried out in the — 100 [/iV] regime, corresponding 
to a Josephson frequency of — lO^'^ [Hz]. 

We now have for the tunnelling Hamiltonian and cur- 
rent operator in the interaction picture 



HT{t) = T{t)+T\t) 
lB{t) 



iKnt 



Ine 



-iKnt 



(32) 
iQe{T{t)~T^t)) (33) 
T(t)==^r,e^-«*V(x„2/„i) (34) 

i 

V{x,y,t) = VH,{x,y,t) = e'^«*Vs(x,2/)e-*^°* (35) 

The effect of the DC voltage on the time evolution of 
the tunnelling operators V is completely captured by the 
phase factor e"^'^*. 

The effective replacement of the tunnelling coupling 
constant F by a time dependent one, F — >■ Fe*""*, can 
also be obtained by performing a suitable gauge trans- 
formation; one that gauges the scalar potential of both 
edges U to zero^''*"'^^. Since the quasiparticle operators 
ip are charged, the tunnelling operators T pick up a phase 
term e'"Q* under this gauge transformation^^. 



IV. LINEAR RESPONSE OF THE 
TUNNELLING CURRENT 



D. Time evolution due to applied DC voltage and 
gauge invariance 

In our approach a simplification is possible which eluci- 
dates some of the later manipulations. In the interaction 
picture the time dependence of the operators, (28), fol- 
lows from the edge Hamiltonian A'o, which includes the 
effect of the DC voltage bias. Since the charge operators 
Qr/l commute with the Hamiltonian Hq we can further 
factorize the time evolution operator as 



-iKot 



g-iHot^ieULQLt^ieURQRt 



The time evolution of the tunnelling operators V due to 
the applied bias voltage can now be made explicit. We 
use the commutation relations of the charge operators 
(11) and the form of the tunnelling Hamiltonian (18). 
This gives for the tunnelling operator V{x^y), 

X Vff„(a;,2/,t)e''=(^^2^+c/flQ«)t ^ (30) 



In the absence of the tunnelling Hamiltonian the tun- 
nelling current vanishes, so {Ib)o = 0- The linear re- 
sponse (29) for the tunnelling current (21) is therefore 

Ib{loq) = {ism = -i I dt{[iB{0),HTmo ■ 



We plug in the expressions for the tunnelling Hamilto- 
nian (32) and the tunnelling current (21) in terms of the 
tunnelling operators T. This gives 



lB{^Q)^Qe di([T(t),rt(0)])o 



(36) 



The correlators of the type (TT) and (T^'T'^) vanish, as 
they describe overlaps of states with different electric 
charge. Furthermore, time translational invariance al- 
lows us to rewrite {[T^t),T{0)])^ = - {[T{-t),T^O)])^. 
A change of integration variable finally results in (36). 

Next we express Eq. (36) in terms of the tunnelling 
operators V{x,y) by substituting Eq. (34) for T. For 



that we introduce the tunnelling-tunneUing correlators 
between the i'th and j'th point contact 



where we replace a^ 



$, 



an with an the relative 



G>{t) = {Vix,,y,,t)V\x,,y„0))o 



(37) 



This gives 



([r(t),rt(o)])o = 5]r,r*e-«* [g> (i) - G<{t)] 



where F^ is the tunnelling coupling constant of the i'th 
point contact. Inserting this into the expression for the 
tunnelling current, (36), the integration over time results 
in an expression in terms of the Fourier transform of the 
G-correlators 



lBiojQ)^Qe^h,iujQ) (38) 

I,,{ojq) - F,r* [G>{uq) - G<{uq)] . 

A final simplification can be made by making use of com- 
plex conjugation, which relates Gfj{uj) = [G|'j(a;)] , and 
the Kubo-Martin-Schwinger condition^^. The KMS con- 
dition applies to two-point equilibrium correlators and 
relates (i(t)B(0))o = (S(O)i(i -t- i/3))o. When applied 
to the tunnelling-tunnelling correlators G we obtain 



G>{t) = G<{t + iP) 
G>{Lo) = e^-G<{uj), 



T^O 



(39) 



and so 



hjiujQ) + IjiiuQ) = 

2|F,F*|Re[e'"'^(l-e-^"«)G>(c.Q)] . (40) 

Here we introduced a.y as the relative phase between 
the couphng constants F^F* = |FiF*|e*"'J. One contri- 
bution to this phase is the Aharonov-Bohm (AB) effect. 
Quasiparticles traversing along different point contacts 
enclose a different amount of flux, which causes an AB 
interference. This interference is independent of the ap- 
plied voltage, provided the geometry is fixed as a func- 
tion of this DC voltage''*; an assumption which does 
not always apply. We define $q — h/{Qe) as the unit 
flux quantum for a particle with Q. The enclosed flux 
quanta between two point contacts i and j is then given 
by ^ij = 27r($i — $j)/$q, where $i is the total flux en- 
closed by the path of quasiparticle tunnelling along the 
i'th point contact. We have for the tunnelling current 



N 



Ib{uq) = Qe(^ |F,|2 (1 - e-P-) G>{uq) + 

i=l 

N 

2^] |F,F*| Re [e'*..+».. (f - e-^-«) G>{u;q)] ) (41) 



phase of the point contacts. The first summation is the 
sum of the tunnelling current through each point contact 
in the absence of any interference. All interference effects 
are encapsulated in the second summation, which we call 
the interference current. 



V. CORRELATORS 

The tunnelling current is completely determined 
through the G^ correlators. In terms of the quasipar- 
ticle operators (10) these correlators are given by a prod- 
uct of four-point correlators, one correlator for each edge 
channel, 



G>(t) = {^PHy^,t)^p{x,,t)^Hxj,0)^P{y,,0))o 
X {a{yi,t)a{xi,t)a{xj,0)a{yj,0))o 



^i^v(Vi:t)^-t^vix,.t)^t^V{xj,0)^-i^v{y,,0). 



(42) 



A. The neutral mode and conformal blocks 

As it stands, the correlator for the neutral channel as 
stated in Eq. (42) is not uniquely defined. Non-Abelian 
quasiparticles span an internal, non-local Hilbert space. 
This is the realization of the non-Abelian statistical prop- 
erties. In the language of conformal field theory^^'^^ 
this internal space is identified as the space of confor- 
mal blocks and the correlator (42) is a particular vector 
in this space. To identify this vector we first need to 
choose a basis in this space of conformal blocks^^'^"'*"'*^ 

The conformal blocks in the correlator correspond to 
the different, possible fusion channels of the quasiparti- 
cles a and a. Symbolically the fusion rules of the fields 
a and a are indicated as 



axa = Y,N',,t 



(43) 



The sum runs over all primary states 9 or quasiparticle 
types of the corresponding chiral algebra, including the 
vacuum state. The integers N^g > are non-zero when- 
ever a field 9 is present in the fusion channel of a and a. 
This fusion rule signifies the possible outcomes when the 
two quasiparticles, a and ct, are brought in close prox- 
imity. In this limit the quasiparticles fuse together and 
either form a new quasiparticle or they annihilate to the 
vacuum. Generally, a correlator such as Eq. (42) repre- 
sents a superposition of possible fusion outcomes. This 
superposition is determined by the history of the system. 
More concrete, the correlator is a linear combination 
of conformal blocks, where each conformal block corre- 
sponds to an intermediate fusion channel. We write sym- 
bolically 



{craaa) 



agSg 



(44) 



i<j 



The sum runs over those primary fields 9 which appear 
in the fusion channel of a and a. With our choice of 
a and a there is always one channel that corresponds to 
the identity or vacuum channel. The functions Eg are the 
conformal blocks and depend on the coordinates of the 
quasiparticlcs. The coefficients ag do not follow from the 
correlator itself but are determined by the history of the 
quasiparticles. 

This summation already assumes a certain order in 
which the quasiparticles are fused together when the cor- 
relator is evaluated. This order is essentially a choice of 
basis in the space of conformal blocks. A different or- 
der in which the quasiparticles are fused together cor- 
responds to a different basis. The corresponding basis 
transformation that relates the two bases is determined 
by an object known as the i^-matrix^"^ . To compute a 
four point correlator, such as G^, we therefore need to 
choose a suitable basis of the space of conformal blocks 
for which the coefficients ag are known. 

In the case of the G^ correlators the quasiparticles 
are formed from the vacuum in pairs at a point con- 
tact. This means the initial fusion channel is the vac- 
uum channel with respect to this basis. Put differently, 
the tunnelling operator V{xi,yi) creates a quasiparticle- 
anti-quasiparticle pair from the vacuum at the i'th point 
contact. It is therefore natural to use this basis, as the 
correlator is a single conformal block with respect to it, 

{a(yi,i)d{xi,t)a{xj,'d)d{yj,Q)) = f vac • (45) 

Pictorially we have^^'*^ 



a{xi,t) cf{.Xj,Q) 



f — 



^■(y*,*)- 



ct(%,0) 



We now identified the vector in the space of conformal 
blocks corresponding to the G^ correlator. However, a 
problem with this basis is that it makes use of fusing 
quasiparticles on different edges. The conformal block 
5vac has components which corresponds to overlaps be- 
tween the two edges. We need to project out these over- 
laps, before explicitly calculating the correlator*^'^^. 

To perform this projection, we switch to a basis in 
which we first fuse together the quasiparticles on the 
same edge, followed by fusion of the these fusion prod- 
ucts. We have 



(cr(2/i,i)CT(xi,t)cr(a;j,0)cr(yj,0)) 



Note that the quasiparticles of each edge are paired to- 
gether and in particular the vacuum channel is always 
present. The coefficients Ovac and ag follow from the 
basis transformation which relates the blocks T and £, 
and they are determined by the components of the F- 
matrix^^. In particular. 



^F 



a a 
a a 



(47) 



vac, vac 



All conformal blocks Te as appearing in Eq. (46) with a 
fusion channel different from the vacuum (Q ^ vac) van- 
ish in the large system-size limit. This is the limit in 
which the size of each edge is taken to infinity, but where 
the distance between the point contacts is held fixed. The 
conformal block that remains corresponds to the vacuum 
channel, and it factorizes into a product of two-point cor- 
relators. We have £vac — a vac -7>ac + • • • and so 



(cr(2/,,, i)cr(xi, t)cr(a;j, 0)cr(yj, 0))o 

= avac(cr(2/i, ^^■(yi, 0))o(CT(a;j, i)cr(a;j, 0))o 



(48) 



The dots represent finite-size corrections which will be ig- 
nored. The two-point correlators are non-zero only when 
a and a fuse to the identity, which is why we started with 
this assumption. What we have accomplished here is a 
disentangling of the edges. In this basis the projection 
onto well-separated edges can be performed. 



B. Two-point correlator of a conformal field theory 

Two-point correlators in a conformal field theory are 
strongly constrained due to symmetries of the CFT^^'*^. 
Following Ref. 79 we first consider the two-point correla- 
tor of some quasiparticle (primary) operator O, 



(o(zi)a(z2)) = 



1 



(Zi - Z2)3 



(49) 



Here the Zi are complex coordinates of the plane, the 
parameter g is called the algebraic decay and it is related 
to the scaling or conformal dimension h of the field O and 
O through g = ^h. The fields O and O must have the 
same conformal dimension or else the correlator vanishes 
identically. 

A temperature is introduced through the conformal 
mapping of the plane to the cylinder, given by z = 
exp(27TiTw/v) where T is the temperature of the sys- 
tem, V is the velocity of the channel and we work in units 
where fc^ = h = 1. The fields transform covariantly^^'^^ 



c^vac + Y^ aeFe (46) according to 0{w) = (j^)^ 0{z), which leads to 



where the basis is now given by 

^iVj^O) (j{xi,t) 
■^0^ a{y,,t) \ ^ \ a{xj,Q) 



{0[wi)d{w2)) = 



{tiT/vY 



sin(7rT(wi — W2)/v)s 



(50) 



This transformation introduces a compactification of the 
coordinates, which is a geometric realization of the tem- 
perature. The Euclidean-time expression is obtained 



through the relation w = vt ± ix. The sign choice de- 
termines the chirahty of the CFT, and a minus sign (— ) 
results in a right moving channel. The real-time expres- 
sion is obtained by performing a Wick rotation. The 
rotation introduces the infinitesimal regulator'", which 
we call S. We have wi — W2 = S + i{vti2 — X12), where 
ti2 — ti ~ ^2- This results in 



gUrQr on the time evolution of the quasiparticle opera- 



{O{x,t)d{0,Q)) 



(ttT/v) 



s[n{TTT{S + i{t + x/v)))9 



(51) 



This correlator is sometimes referred to as the greater 
Green's function. In the end the propagator is neatly 
summarized as 



{Oix, t)O{0, 0)) = v-<^Pgit - x/v) 
1 






sm{TrT{S + it))9 



r = o 



r>o 



(52) 



For completeness, we have included the zero-temperature 
limit. Putting everything together we obtain for the cor- 
relator of the neutral mode 

{a{y„t)a{x,, t)cr{xj,0)a{yj,0))o = 

aviic{cr{yi,t)a{yj,Q))o{a{xi,t)a{xj,0))o^ « 

^Pg„ {t + iia/v^)Pg^ {t - rjb/vn) (53) 

Here a — \yi — yj\, b = \xj — Xi\, and u„ and g„ = 
^hrj are the velocity and algebraic decay of the neutral 
channel. The parameter ry = ±1 denotes the chirality of 
the neutral channel relative to the charged mode, with 
(77 = +) representing the same chirality. 



C. Charged mode 

The charged mode is Abelian, meaning all fusion chan- 
nels are unique and the projection onto disentangled 
edges can be done without having to perform a change of 
basis in the space of conformal blocks. The factorization 
is 

/ i-^univi.i) —i-^-u>(xi,i) i-^iplxjfl) — i-2=(p(u ■ 0)\ 

X (g-^^'^(--*)g^ J;'^(--o))o + . . . (54) 

The dots represent finite-size correction which we ignore. 
It is tempting to apply the results for the two-point cor- 
relator as motivated in Section VB. Although this leads 
to the correct result, it glosses over the fact that the cor- 
relator for the charged mode is taken with respect Kq 
instead of the usual conformal Hamiltonian if g ■ 

At this stage we recall that we have already taken 
into account the effect of the coupling terms —cUlQl — 



tors e 



',f(Vi) 



This was done in Section HID and leads 



to the phase factor g"^Q*. However the coupling terms 
also appear in the initial density matrix, e~' °, which 
can potentially give rise to extra contributions. We now 
show that these contributions are attributed to finite-size 
correction which we ignore. We do this by explicitly cal- 
culating the propagator with respect to Kq using a mode 
expansion. 

We assume a finite system length L with periodic 
boundary conditions''^'^'' on dip and switch to a Fourier 
decomposition for d^ip, 

\j^dMx) = Po-iJ2 ^(^'""h - e-^'^-bl) (55) 



fc>0 



Y ^'/'(a;) = ^o+xpo-Y, TK^""'^^^ 



kxifi ^ 



k+e b^) . 



Here k — ^ixnjL with n integer > 0. The zero mode po is 
proportional to charge operator p^ — ^J2ti /vLQ. From 
the commutation relations (5) we obtain for the modes 



[•^o,Po] =« 



[&lM = i 



(56) 



and the remaining commutation relations vanish. Up to 
a constant term the normal ordered Hamiltonian is given 
by 



K. 






'^ V 9~'^° + "<= X! k^\^k 



(57) 



fe>0 



From the Hamiltonian we derive the time evolution of 
ip(x,t) with respect to Kq. This gives 

<f{x,t) = -^eUt+ J -j-fo + y -j-{x + Vct)po 
- \[^y ^(e^M-+-ct)5 ^ g-^fe(x+..t)^t ) (58) 

The Hilbert space is constructed in the usual way^", 
meaning we have a vacuum state |0) which satisfies 
Polo) = bk\0) = 0. The modes bl with fc < act as 
creation operators of momentum modes on this state. 
Since ipo does not enter the Hamiltonian, the opera- 
tor po is conserved and can be diagonalized simultane- 
ously with the Hamiltonian. The operator e^°"^° creates 
a charged eigenstate, as it raises the eigenvalue of po by 
a, i.e. po{e"^'^°\0)) = a(e'"'^°|0)). In the large system- 
size limit the overlap between states with different charge 
a vanishes^^. The allowed values for a which construct a 
state with non-zero norm depends on the chiral algebra. 
The normal ordered exponential operator is defined as 



: e 



fc>0 



^-ioiy/veUt^ia,^ ^ipo ^ia^^ ^(x+v^t)pQ 



^fcl-g-^^A 



(59) 
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with the time evolution again due to Kq. We are inter- 
ested in the two-point correlator, 



1 



-Tr 



g-^ArOgiaip(a;,t)g-iaip(0,0) 



(60) 



Note again that the time evolution is with respect to 
Kq. The correlator is computed for each mode separately, 
since the different modes commute and so Z = Zq Ofe ^k- 
The contributions of the non-zero modes to this correla- 
tor is the same as in the zero-bias case, see e.g. Ref. 70 for 
details on this computation or Ref. 83 for an alternative 
derivation. We have 



Tr [e"'^"=''fc'"' e"''"\/?^' 



ifc(x + «et)!,t 






'PgSt + x/Vc) + .. 



(61) 



The normalization of the two-point correlator is unity 
in this limit. For the zero mode we first note that com- 
pactification of the boson ip = if + 2tiv restricts the spec- 
trum (eigenvalues) of po to J^ with k integer. This 
gives for the correlator, 

1 



^f-l^-Pi^ Pa- V^eU ^J j^ po) i^-iaV^eUt ^ia^ '^ifo 



iay/^(x+v^t)pa -ia^- 



— ioi^/ueUt I 



The dots represent finite-size corrections. The effect of 
the zero mode coupling to the external potential on the 
two-point correlator is the phase factor e-'"v^eC/t_ -^g 
obtain for the two-point correlator with the time evolu- 
tion due to Kq 



1 



-Tr 



Q-l3Ko^iaip{x,t)^- 



11^(0,0) 



^-icx^eUt^- 



"^Pg^{t + x/vc) . (62) 



Finally, for the expression for the desired correlator 
appearing in Eq. (54) we strip off the phase factor 
^-la^eut ^ since we already extracted this through the 
manipulation performed in Section III D - it leads to the 
phase factor e"^Q*, which is already taken into considera- 
tion in the expression for the tunnelling current (41). We 
obtain 



/g*;%'^(!'»'*)g^*;%'^(^»'*)g*;%'^(^3'0)g~*;%'''(w^o)\„ — 



'v^'' 



;o = 



v-^sap^^ (i _,_ a/vc)Pg^ (t - b/vc) + 



(63) 



This form matches with what we obtain by simply replac- 
ing the two-point correlators (54) by the propagators Pg. 



D. Quasiparticle braiding and bulk-edge coupling 

The correlators of the neutral and charged modes, 
equations (53) and (63), encapture part of the dynamical 



effects of quasiparticles traversing along the edge. The 
other dynamical contribution is due to the AB phase. 
In addition, there is also a topological contribution to 
the tunnelling current due to braiding of bulk and edge 
quasiparticles'*^'^^'^*"^^. The correlator G^ is interpreted 
as the amplitude of the process in which a pair of quasi- 
particles ip and "0^ are created from the vacuum at the 
j'th point contact and annihilate to the vacuum at the 
i'th point contact. If one or multiple quasiparticles is 
present between these point contacts, the resulting am- 
plitude contains a contribution coming from the quasi- 
particle braiding. This so-called matrix element is de- 
picted in Figure 3. 

More generally. Figure 3 represents the expectation 
value of Wilson lines computed with respect to the full 
topological quantum field theory and it is fully deter- 
mined in terms of the 5-matrix^^. To fully determine this 
expectation value we require to specify the exact TQFT 
and the configuration and state of the bulk quasiparticles. 
In general the outcome is some complex valued function 
Aij{x): bounded by |-4ij(x)| < 1, which depends on the 
topological quantum number x associated with the bulk 
anyons inside the interferometer. For the G^ correlators 
we have 



Gfj — Ovac Aij (x) X (dynamical contributions) 



finite-size effects . (64) 



The effect of quasiparticle braiding is a topological ef- 
fect, due to the statistical properties of the anyons. In 
the case of the Moore-Read state the effect leads to what 
is known as the even-odd effect'*^'^^'"'®"^^. When there 
are bulk quasiparticles present inside the interferometer 
and these quasiparticles are located far from the edge 
then the interference current due to tunnelling of the e/4 
quasiparticle vanishes when the number of bulk quasipar- 
ticles is odd. When the number is even the interference 
current re-emerges. 

The situation is more complicated when the bulk quasi- 
particles are close enough to the edge of the system. 
In that case the coupling between the bulk quasipar- 
ticles and edge degrees of freedom needs to be taken 
into account^^"^^. This coupling can induce tunnelling 
of the neutral degrees of freedom associated with the 
non-Abelian statistics from the bulk quasiparticles to 
the edge theory. One result is that even in the case 
of an even number of bulk quasiparticles located inside 
the interferometer this bulk-edge coupling can effectively 
flush out the interference current. Averaged over time 
the tunnelling of neutral degrees of freedom can greatly 
reduce the strength of the interference current. We do 
not take into account the effect of bulk-edge coupling, 
but we do note that this effect can be relevant to recent 
experiments'^ '^''^^''^^ 
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Figure 3. Quasiparticles inside the interferometer braid with 
quasiparticles tunnelling along the point contacts. At lowest 
order the effect of braiding is captured by the corresponding 
braiding diagram, which is determined from the topological 
quantum field theory. 



E. G^ correlators and its Fourier transform 

The expression for the G^ correlator (42) follows 
straightforwardly from combining the correlators for the 
neutral and charged mode, (53) and (63). 

avac(V'^(yi,*)^(2/j,O))('0(xi,t)V''^(a;j,O)) + . . . 



om 



PsAt^^JPsAt-v^) 



(65) 



Xj I as the 



Here we have defined a = \yi — yj\ and b = \x 
distance between the i'th and j'th point contact along 
the upper and lower edge respectively. Recall further- 
more that ?7 = ± represents the chirality of the neutral 
channel relative to the charged channel and Aij (x) is due 
to braiding of quasiparticles. 

For the tunnelling current we need the Fourier trans- 
form of the G^ correlator. In Appendix A we show how 
this Fourier transform is obtained. We first treat the con- 
tribution due to tunnelling along a single point contact, 
Gf^ . The correlator for Gf^ (t) is independent of position, 
since a = & = in (65). We have 

G>it) = v~'a-v~^s^P2,{t). (66) 

with g — gn + gc twice the total scaling dimension of the 
quasiparticle. Using the result of (A3) gives 

G>i-)=e^^^^^B(g + ^A^,g-^A^) (67) 



23„ 2ge 
Vn Vc 



27rr'' 



27rr/ 



Here B{x,y) is the Euler beta function and we have set 
the integral regulator S to zero. We treat the zero tem- 
perature case later on. 

The expression for the more general case {i ^ j) is 
more complicated. We write the Fourier transform of 
GfAuj) {i ^ j) as the integral definition of Carlson's R 
function"''^. This is a multivariable generalization of the 
Gauss hypergeometric function. An alternative way of 
representing Carlson's R function is through the fourth 
Lauricella hypergeometric function^^"^^ , see also the ap- 
pendix. We have cf. Eq. (A6) the following expression, 

(68) 



GfA^Q\ 



^ am)h7;°\^q)gi{ujq) 



where all (trajectory-dependent) interference effects are 
hidden away in the modulating function H™°'^ given by 



h: 



mod / \ 

._OJQ_ 



xi?(ff 



^ ^^T(b-a)(fj+,,2^) 
;{5c,5'c,5'n,5'n}; 



-2ttT^ 2TrT-!^ -r727rT-5- r72irT-!^ 

e "c , e "<: , e ' "" , e ' 



(69) 



The R function is treated extensively in Ref. 57 and we 
have summarized some of its properties in Appendix B. 
In particular, the order in which the parameters appear 
in (69) is relevant for its evaluation. Furthermore the 
R function allows for certain transformations of the ar- 
guments, see also the appendix. Computation of the R 
function is explained in Appendix B 3 using results of 
Ref. 84. We mention one transformation in particular 
which is equation (B5). Through this transformation we 
have the equivalent expression of the modulating func- 
tion (69). This transformation effectively switches a O 6 
in the expression of Hf^°'^ and simultaneously changes 
the sign of ojq , 



xPi^g + i^^'Agc, gc, gn, gn}; 



2-kT^ -2irT-!^ r727rT^5_ ~n2TTT-^ 

e "c , e "= , e ' "" , e ' 



(70) 



The function Aij (x) describes the effect of possible quasi- 
particle braiding entering the correlator Gij. Finally, in 
the expression for GfJuj) we recover the expression for 
the single point contact case, Eq. (67). The effect of 
the spatial separation of the point contacts, and thus all 
interference effects, is completely captured by the mod- 
ulating function i/™°'^. 

Since the R function is so closely related to the Lau- 
ricella function we also mention the form of the Fourier 
transform in terms of this function. The exact relation is 
explained in the appendix. Here we assume for simplic- 
ity a symmetric interferometer b = a. Assuming -5^ > -5l 
and the expression reduces to 



rrniod 



(wq) 



— 2nT-^-n iuJn- 



.. 77(3)/ ■ '^Q f 1 o 1 -2TrTa(J- + i) 



^_g-2.Ta(i-i)^^ 



-4ttT- 



(71) 



This expression no longer depends on the chirality pa- 
rameter r] — ±. The symmetric interferometer does not 
distinguish between chiral and anti-chiral edge states. 



VI. EXPRESSION FOR THE TUNNELLING 
CURRENT 

We combine the expression for the tunnelling current 
(41) with the expression for the Fourier transform of the 
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correlators (67) and (68) and obtain the following expres- 
sion 



Ib (wq) = 2Qe^2^^-^^avac|rcff(tJQ)|' 



(^\ 



\2TJ V 27rT'^ 2ttTJ 



) ■ (72) 



In the spirit of Ref. 40 we have combined the effects due 
to interference into an effective tunnelling coupling am- 
plitude, 



AT 



i=l 
N 

2^ |r,r,| Re[A,(x)e^*-+'"--i?,7'^(a.Q)] . (73) 

The function H^°'^{luq) is given by (69), which we call 
the modulating function. The r,;'s are the tunnelling 
coupling constant of the i'th point contact and Uij is the 
relative phase between F^ and Tj. We also introduced 
the Aharonov-Bohm phase ^ij and contributions due to 
quasiparticle braiding are attributed to Aij (%) . The dis- 
entangling of the conformal blocks results in the factor 
flvac- Within our setup only i?™°'^ depends explicitly on 
the external voltage bias. The tunnelling constants Ti de- 
pend on the exact geometry of the interfcrometric device, 
and so the normalization of the current is not universal. 
Expression (73) for the tunnelling current is of the 
form, 

N 

/i3(^Q)- (E|rd'+E^r')^^inglcpc.(^) (74) 

i i<j 

= Iq + lose 

F^°'' = 2|r,r,| Re[^,,(x)e'*'^+^"-i/,7'^(^Q)] 

All interference effects are contained in the function 
iT'mod^ which we call the interference term. We only deal 
with interference between pairs of point contacts; there 
are no interference effects involving tunnelling along 
three or more point contacts. This is due to the linear 
response approximation, which only takes into account 
effects up to order Ir^rjl. 

The modulating function Hf^°'^ is a function of the 
different energy scales, which are set by the temperature 
and voltage bias, and the scales associated with the ve- 
locity and distance between the point contacts. 
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These parameters enter the expression for the function 
H^°'^ through dimensionless combinations, and the func- 
tion depends on the relative scales. The modulating func- 
tion is, up to an exponential factor, determined by Carl- 
son's R function which we treat in the appendix. The 



R function is a scaling function, which manifests itself 
through the homogeneous scaling transformation (B4). 
It is computed through its relation to the Lauricella func- 
tion and the corresponding Taylor series as described in 
Appendix B 3. 

The expression for the interference current is very gen- 
eral, and the price we pay for this is a limited intuition 
when it comes to the behaviour of the corresponding 
modulating function, H^°'^. We can still summarize the 
general behaviour of the function as a function of the 
physical parameters. As a function of increasing volt- 
age ojQ the modulating function is the sum of multiple, 
decaying oscillations. The frequencies of the oscillations 
are determined by the edge lengths and edge velocities. 
The temperature and algebraic decay determines the rel- 
ative amplitudes of the oscillations. In addition, for large 
temperatures _ff ™°'^ decays exponentially. Some of these 
features are proven analytically, while others follow em- 
pirically from numerical analyses. 



VII. SPECIAL CASES AND 
GENERALIZATIONS 

The main result of our work is the expression for the 
interference term (73) for the tunnelling current (72) in 
terms of the R function (69). Here we consider several 
limits and generalizations, such as the zero temperature 
limit and other cases in which the expression for the mod- 
ulating function H™""^ simplifies. This relates our results 
to earlier work^"'"'^"^'''''^ . We consider the generalization 
to more than two modes and discuss a relation to the 
two-point quasiparticle propagator. 

Recall that we use gc and g„ to denote the algebraic 
decay of the charged and neutral channel, and v^. and 
Vn the corresponding edge velocities and ry = ± as the 
chirality of the neutral mode. In the case of three or 
more point contacts we obtain a modulating function for 
each unique pair of point contacts, Hf-°'^. We use a and 
6 to denote the length between the I'th and j'th point 
contact along the upper and lower edge respectively. In 
principle, these lengths depends on i and j, so a = aij 
and h — bij. However, we omit these subscripts for the 
sake of breviety. 

Finally, we set g — gc + g-n as the total algebraic decay 
and work in units where ks — ti ^ 1. 



A. Zero temperature limit 

The zero temperature limit can be obtained in two 
ways. The first is to start with the expression for the 
propagator at zero temperature, (52), and follow the 
same steps as in the finite-temperature case by comput- 
ing the Fourier transform of G^ and G"^. Alternatively, 
we can start with the expression for the tunnelling cur- 
rent at finite temperature, and from here take the zero 
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temperature limit. Both routes should produces the same 
result. 

However, the first route leads to an obstruction. When 
we attempt to determine the interference current we en- 
counter the following integral (see also Appendix C) 



GfJu;) ~ G<iuj) 
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We do not know how to solve this integral with these 
general parameters and we are not aware of a reference 
in which it is treated. Therefore we proceed with the 
other route, in which we start with the finite temperature 
expression, Eq. (72), and take the zero temperature limit. 
For the current we find the usual power-law behaviour 
times an effective coupling amplitude 



Ib (wq) 



27r 



2Qe avac|reff(c^Q)na;Qpg-^sgn(c.Q) . (77) 

Vn Vc 

The expression for |rcff(LL'Q)p is the same as in the finite 
temperature case, Eq. (73), but with a different expres- 
sion for the modulating function H^°'^. We have worked 
out the zero temperature limit of H^°'^ in Appendix C. 
The result is 
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The function $2 is the confluent Lauricella hypergeo- 
metric function of 3 variables^^ and its series representa- 
tion is given by Eq. (C18). It can be extended to include 
more than two modes per edge. This expression for Hf^°'^ 
should also be obtained by direct computation of the in- 
tegral (76). 

In the symmetric case where a = b the modulating 
function reduces to 
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The function $2 is known as a Humbert confiuent hy- 
pergeometric function of two variables^'^'*^. For the sym- 
metric interferometer the chirality of the neutral mode 
has no effect on the current. 



effectively an edge with a single channel. The modulating 
function Hl^°'^{uj) reduces to the Gauss hypergeometric 
function. 
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The function 2F1 is the Gauss hypergeometric function. 
For equal distances between the tunnelling points on both 
edges (a = b) the expression coincides with that found in 
Ref. 40, although to arrive at this expression we require 
some manipulations of the Gauss function. These can be 
found at e.g. [p. 1009] in Ref. 60. 



eq. (79) = 27r 
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Expression (80) seems obscure and overly complicated in 
comparison with (79). However, the representation (80) 
is an expansion in terms of the parameter e"'^'"^^' , which 
tends to zero for large temperature. In contrast, the ex- 
pansion appearing in (79) is in terms of 1 — e~'^'^'^^^^ , 
meaning the argument of the Gauss function tends to one 
for high temperatures. The exact behaviour of the Gauss 
function around unit argument is problematic, and leads 
to slow convergence of its Taylor series or even singular 
behaviour. In fact, the standard way of analysing the 
behaviour of 2-^1 (a, b;c;l — z) for z ^- is by first trans- 
forming it into a function of the form 2-^1 (a', &'; c'; z). 

The zero temperature limit can again be obtained in 
two ways: by directly computing the Fourier transform or 
by taking the zero temperature limit of the finite temper- 
ature expression. In this case it is possible to determine 
the Fourier transform directly, which we have done in 
Appendix C 1 . We also show that this Fourier transform 
matches with the zero-temperature limit, demonstrating 
the equivalence of both routes. We find 
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Here J„_i(x) is the Bessel function of the first kind and 
y 2 

this expression matches with what was found in Ref. 40 
when we set a — b. 



B. Equal velocities and chiralities 



C. Fast charged channel 



For equal velocities and equal chiralities between the 
two channels we set v = Vn = Vc and rj = +1. This is 



We consider the limit where the energy scales asso- 
ciated with the charged mode are far greater than the 
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remaining energy scales, 
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The scales on the right hand side are that of the neutral 
mode, the temperature scale and the applied voltage bias. 
In this limit the modulating function is 
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On the final line we obtain the first Appell hypergeomet- 
ric function of two variables^*^ , Fi(a; /3,7; zi, Z2). When 
f/c = this function reduces to the case of a single edge 
mode Eq. (79), as expected. 



D. Large interferometer and high temperature 
limit 



For well separated contacts we consider large a + b. 
In Appendix B 2 we show how this behaviour can be 
extracted from the integral. This limit suppresses the 
interference current exponentially according to 
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This is interpreted as an effective dephasing length 
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Beyond this scale the interference current is suppressed 
as / oc g-i°'+i')/LT -^^itii a + b the total circumference of 
the interferometer. A similar analysis applies for high 
temperatures. Setting 
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and the interference signal vanishes as / ex e^^'^^ . In 
general the decoherence effects are reduced by decreasing 
the temperature. See Ref. 44 for further discussion on 
energy scales and visibility of the interference signal. 



E. Asymmetric interferometer 

We now consider the limit where the length of one edge 
approaches zero. We set a = which effectively merges 



the point contacts on one edge. We obtain 
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The reduction of this expression to the corresponding hy- 
pergeometric form depends on the sign of 77. For 77 — +1 
we have 









while for 77 = —1 we obtain 
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Here Fi is the Appell hypergeometric function of two 
variables^°. Using transformation properties of the Ap- 
pell function, which can be found in e.g. [p. 1020] in 
Ref. 60, for the case of 77 = —1 we can obtain a single 
expression for the H^°'^{uj) function given by 

T-, ( ■ ^ r, 1 -27rT^ 1 -27rTr7^!-\ 



F. More than two channels 

Our result for the interference current generalizes to 
edges which consists of more than one mode, all with dif- 
ferent velocities. We can also include the possibility of 
different edge velocities for each edge. The edge veloc- 
ity is not a topologically protected property of the edge 
mode, and its value(s) can depend on the exact geometric 
details of the corresponding device. 

The generalized result is obtained if we assume the 
modes decouple in a similar fashion as in the two-channel 
case or that the propagator factorizes along the lines of 
(65). In these cases the correlator G^ generalizes to 

ra 

G> (t) = avac^»j (x) n \l\r Pg, (^ + majv^^L) 

xPg^{t-iq,b/v,^R) (89) 

Here gi and 77^ are the algebraic decay and chirality of the 
j'th edge channel and Vi^ji and Vij^ the velocity of the i'th 
edge mode on the lower and upper edge. The function 
Aij{x) accounts for possible braiding of quasiparticles 
and Ovac arises due to disentangling of the edges. The 
current is still determined by the Fourier transform, and 
the only change arises in the modulating function and 
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the normalization of the tunnelhng current which now 
involves all of the velocities, see (72). We have 
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Here the arguments are ordered sets consisting of the 
algebraic decay and energy scales, 

{9i,9i}7Ll = {91^91,92,92,- ■■ ,9ni,9m} 
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and g — X^iff*- Computation of this function is similar 
to the two-channel case and covered in Appendix B 3. 



G. Two-point correlators and the R function 



The tunnelling correlator Gf^ is constructed through 
projection onto decoupled edges, which results in a de- 
composition in terms of a product of two-point correla- 
tors. A simpler expression arises when we consider the 
the two-point propagator of a non-Abelian anyon on a 
single edge. We have in the conformal limit 



X>^^^{t,x) = {^P\x,t)^iO,0)) 
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Here we absorb the chirality of each mode into the veloc- 
ity Vi, which can therefore take on negative values. The 
corresponding Fourier transform with respect to time is 
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where g = J2i 9i ^-^d all spatial dependence is captured 
by the function 



H(„i){uj,x) =e 
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This is the equilibrium two-point quasiparticle propaga- 
tor in a frequency-coordinate representation. 



VIII. PLOTS OF THE MODULATING 
FUNCTION AND INTERFERENCE CURRENT 

In this section we plot the modulating function 
and the corresponding interference current. Based on 



experiments'^"'^'^^''^'*"'^^ we take the distance between two 
point contacts to be around 2 [/xm]. For the velocity no 
experimental data is available, but numerics^® suggests 
a much faster velocity for the charged mode compared 
to the neutral mode on the order of Vc/vn ~ 10 and 
Vc ^ 10"* [m/s]. The applied voltage bias lies typically in 
the range of 10 to 50 [/iV] and temperature ranges in the 
order of 10 — 25 [niK]. We assume a lower temperature of 
~ 1 [mK] as this significantly improves the rate of con- 
vergence of the series used to compute the expression for 
the tunnelling current, see Appendix B 3. 

In this section we are mainly interested in the be- 
haviour of the R function. The factor ytij(x)e'*'3+*"'J 
is due to quasiparticle braiding, the AB phase and the 
relative phase between the tunnelling amplitudes of the 
point contacts. They are assumed to be independent of 
the applied voltage bias and we set the total factor to 
unity. We comment on the AB effect in the next section. 

The final parameters that need to be fixed are model- 
dependent, and correspond to the filling fraction i/, the 
algebraic decay of the quasiparticle propagators 5„ and 
5cj and the quasiparticle charge Qe. For a given edge 
state a renormalization group analysis predicts the quasi- 
particle with the lowest algebraic decay, gn+gc, to be the 
most relevant perturbation'^''''*^'^^. Quasiparticles with a 
larger algebraic decay are less relevant in the language of 
the renormalization group and we ignore their contribu- 
tions in the plots. 

A second effect is that the effective magnetic length, 
1% = h/{QeB), is larger for quasiparticles with a smaller 
charge. The bare tunnelling matrix element depends on 
this length scale, and it is expected that a smaller charge 
correspond to larger matrix elements. Some trial states 
predict multiple quasiparticles with the same algebraic 
decay. In these cases the contributions to the tunnelling 
current is expected to arise from the quasiparticles with 
the smaller charge. 

Computation of Carlson's R function is not completely 
straightforward. The function admits a multivariable 
Taylor expansion or one can resort to numerical integra- 
tion of the Fourier transform G^. Using combinatoric 
results of Ref. 84 the Taylor expansion is cast into a sin- 
gle summation, which we explain in Appendix B 3. We 
use this expansion for computing the R function. 

For physically relevant values of the input parame- 
ters both the series expansion and numerical integration 
schemes converge very slowly. In particular a higher tem- 
perature scale reduces the convergence rate significantly. 
We apply a series acceleration using the CNCT method** 
to partially remedy this problem, see also the appendix. 
However, even the CNCT method is not practical for high 
temperatures and to our knowledge an efficient numerical 
scheme is still lacking. 

Due to these convergence problems we are not able 
to compute the R function for all ranges of the physical 
parameters. For instance, we mostly assume tempera- 
tures of 1 or even [mK] for the sake of convergence of 
the modulating function. We also plot the modulating 
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function over a range of the source-drain voltage which 
hes outside of what is reached in experiments. We have 
chosen for this range as we want to demonstrate the non- 
trivial behaviour of the R function over a greater voltage 
range. Finally, we point out that currently the edge ve- 
locities have not been measured and it is possible that 
the values used in the plots are inaccurate. 



A. The tunnelling current without interference 

Before we provide some plots of the modulating func- 
tion and the interference current we first discuss the gen- 
eral behaviour of the tunnelling current in the absence 
of interference'^^. This expression also enters the result 
for the total interference current. It is given by Eq. (72) 
in terms of the Euler beta and the hyperbolic sine func- 
tion, with the tunnelling amplitude held constant. The 
tunnelling current is characterized by the total algebraic 
decay, and we discuss two particular values of g for which 
the function simplifies. These follow from the properties 
of the gamma function^^ 
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In the limit of ojq — )■ oo and g = h the expression for 
the tunnelling current approaches a constant value, while 
g — 1 grows linearly with ujq. For the remaining cases 
the current decays to zero for g < i^, grows sublinearly 
for ^ < g < 1, and grows superlinearly for g > 1. Finally, 
at zero temperature the expression for the current follows 
the power law behaviour 



/b oc IcjqI^s sgn(cjQ) 



(95) 



while for high temperatures the function follows 

Ib cx u;QT^a-2 _ (gg) 

B. The tunneUing current with interference 

The upper panel of figure 4 is a plot of the total tun- 
nelling current (eq. (72)) with and without interference 
for the case of the Moore-Read / PfafRan quantum Hall 
trial state^^'^^ for the v = 5/2 plateau. The lower panel 
of figure 4 is a plot of the corresponding modulating func- 
tion Yle[Hf^°'^], given by equations (73) and (G9). The 
parameters for the set gc = gn = 1/8 and Qe = e/4. 
This result is also analysed in Ref. 42. See the figure 
caption for the exact values of all parameters. 

The normalization of the current, which is the pref- 
actor appearing in expression (72), contains the tun- 
nelling coupling constants F^. These factors are non- 
universal, meaning the normalization of the current is 
non-universal as well. In Figure 4 the current with- 
out interference is normalized by its maximum value. 
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Figure 4. The tunnelling current through an interferometer 
with and without the modulating function. The current is 
normalized by the maximal value of the tunnelling current 
without interference (max(Jfl(i/™°'^ = 0))). The quasipar- 
ticle is the Qe = e/4 QP of the Pfaffian state which has 
Qn ~ Oc = |. The remaining parameters are Vc = 7 ■ 10^ 
[m/s], «„ = 1 ■ 10^ [m/s], T = 1 [mK], a = 2.0 [/im] and 
b — 1.8 [/im]. All coupling constants are equal. Interference 
effects due to braiding with bulk quasiparticles is absent, i.e. 



(94) weset Aj(x) = 1- 



/^(no interference)/max(/B(no interference)). The nor- 
malization of the current with interference is chosen such 
that when the two currents cross in Figure 4 the modu- 
lating function vanishes H™°'^ = 0. 



C. Voltage and geometry dependent oscillations 
and corresponding frequencies 

The modulating function lie[Hfj°'^] shows multiple os- 
cillations and decays when ]/ — > oo, see the lower panel 
of figure 4 and the upper panel of figure 5. A numerical 
analysis (figure 5) shows that for an asymmetric inter- 
ferometer (a 7^ b) and two different edge velocities the 
modulating function consists of four oscillating signals 
with frequencies 



JXj ,Vi 



Vih 



(97) 



where Xj = a,b and Vj = Vc, Vn- These frequencies can be 
extracted from the (x, t) representation of the tunnelling- 
tunnelling correlators GfJt), see Eq. (65). The peak val- 
ues appearing in this correlator correspond to the fre- 
quencies (97). We also find that the frequencies are inde- 
pendent of the temperature and algebraic decay - these 
parameters only influence the total and relative ampli- 
tudes of the oscillations. In the limit of a symmetric in- 
terferometer (a ss b) the number of contributing oscillat- 
ing frequencies drops from four to two, since fa,vi ~ fb,vi ■ 
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Figure 5. Upper panel: the modulating function for a ficti- 
tious quasiparticle with gc — 1/8 and gn = 1/6 at T = [K] as 
a function of the voltage. Furthermore, Qe = e/4, Wc = 3- 10^ 
[m/s], Vn = 9-10^ [m/s], a, — 4.0 [iJ-ra] and b = 2.5 [/im]. Lower 
panel: the corresponding (windowed) Fourier transform. The 
vertical lines represent the predicted frequency components 
given by QeXi/{vih). 



Figure 6. Upper panel: the modulating function for a quasi- 
particle with gc = gn = 1/8 at T = [K] while varying the 
distance of one edge (0 < a < 3.0 [mih]). The applied volt- 
age is kept constant at 1/ = 60 [/.tV]. Furthermore Qe — e/4, 
Vc = 1 ■ 10^ [m/s], v„ — Q- 10^ [m/s] and b = 2.5 [/irn]. Lower 
panel: the corresponding windowed Fourier transform with 
peaks at the predicted frequencies QeV/{vch) and QeV/{v„h). 



In this regime the two oscillations form a modulating sig- 
nal with 'fast' and 'slow' frequencies 2F^(i7" ^ ^)' which 
was also found in Ref. 42. It is also possible that the edge 
velocity for each channel is different on opposite edges. 
In that case we still have four different frequencies in the 
Fourier spectrum, even in the case of a symmetric inter- 
ferometer. 

The second analysis we perform looks at the oscillat- 
ing behavior of the modulating function as a function of 
the length of one edge, while keeping all other parame- 
ters fixed. These are the oscillations in H"^°'^ when a is 
varied. The frequencies of these oscillations are obtained 
through a numerical Fourier transform, see figure 6. The 
modulating function shows a similar decaying, oscillat- 
ing behaviour as in the case of varying the voltage, with 
frequencies given by 
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Since the other edge length b is kept constant we observe 
only two contributing frequencies. For the case of a single 
edge velocity these frequencies can be extracted from the 
expression of the current at zero temperature, (81), by 
making use of properties of the Bessel function. However, 
we are not able to extract the frequencies in expression 
(98) analytically for the more general case. We suspect 
that such a result can be obtained from the R function 
through an asymptotic expansion, which we leave as an 
open problem. We expect that these results carry over 
to the more general case of several edge channels and 
different velocities, see Section VII F. 



D. Effect of temperature 

In figure 7 the modulating function is plotted for the 
temperatures [mK], 10 [mK] and 18 [mK]. The T = 
case is computed using the confluent Lauricella hyperge- 
ometric function, as explained in Appendix C. Compu- 
tation of the confluent Lauricella function is very similar 
to the finite temperature case. 

The convergence of the series representation used to 
compute R function becomes progressively worse for tem- 
perature scales larger compared to the remaining energy 
scales. Computing the R function using the series expan- 
sion in this regime becomes impracticable, even when we 
employ a series acceleration. This type of slow conver- 
gence is similar to that exhibited by the Gauss function 
2Fi{a,b;z) when \z\ — )■ 1. For the Gauss function a set 
of linear transformations exist which allow one to avoid 
this |2:| = 1 singularity^^, see also Eq. (80) and the cor- 
responding discussion. We are not aware of a general- 
ized type of transformations applicable to the R func- 
tion. Due to this slow converge for high temperatures 
we frequently put T = or T = 1 [mK] throughout this 
work. 

From figure 7 we observe that the oscillations are in- 
dependent of the temperature. Other numerical analyses 
suggest that this remains valid for other physical param- 
eters as well. Instead the temperature appears to be re- 
sponsible for the relative and absolute amplitudes of the 
oscillations which were studied in the previous section. 
In particular, higher temperatures cause an exponential 
suppression of the function as was found in Section VII D. 
Lower temperature increase the visibility of the interfer- 
ence signal. 
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Figure 7. The modulating function at three different tem- 
peratures, keeping all other parameters fixed. The tunnelling 
quasiparticle is the e/4 quasiparticle of the PfafEan state with 
9c~ Qn ~ \- The remaining parameters are Vc = 9-10'' [m/s] , 
Un = 9 • 10'^ [m/s], a — 3.5 [/im] and h = 3.5 [/im]. 



in all cases of a decoupled neutral and charged chan- 
nel as described in Section II. In the case of the Anti- 
Pfaffian the neutral and charged channels have opposite 
chiralities. All of these edge theories predict a charge 
of Qe — e/4 associated with the quasiparticle with the 
lowest algebraic decay. Furthermore, the quasiparticle 
with second-smallest algebraic decay is for all cases a 
Laughlin-type anyon with a charge of e* = e/2 and alge- 
braic decays of gc — 1/2 and (/.„ = 0. Figure 8 is a plot of 
the corresponding modulating functions for the different 
edge theories, including the e/2 quasiparticles. 

As we mentioned before, in the language of the renor- 
malization group the most relevant tunnelling operator 
correspond to quasiparticles with the lowest algebraic de- 
cay. In the case of the Anti-Pfaffian the lowest algebraic 
decay is given by \ and it corresponds to two quasipar- 
ticles, the e/4 and e/ anyon. In this case we also need 
to take into account that quasiparticles with a smaller 
charge have a larger magnetic length, and therefore a 
larger bare tunnelling amplitude. So also in the case of 
the Anti-Pfaffian it is expected that the interference cur- 
rent is due to tunnelling of the e/4 quasiparticle. 



The experiments are typically performed at tempera- 
tures of T = 25 [mK] or lower. Numerically, we have 
not been able to reach temperatures higher than T = 20 
[mK]. We expect that the behaviour of the modulating 
function as predicted by our results remains valid in this 
regime. In particular we expect that the frequencies of 
the oscillations (97) and (98) are independent of the tem- 
perature, although we have not proven this analytically. 
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Table I. Parameters for quasiparticles for different edge mod- 
els of the If = \ state. Listed are the quasiparticle charge 
e* = Qe, the algebraic decay of the quasiparticle's neutral g^ 
and charged gc channel, the total algebraic decay g ~ gc + gn 
and the chirality. The e/2 Laughlin quasiparticle is present 
in all three states. Data obtained from Ref. 44. 

The most prominent state for which the correspond- 
ing topological phase is conjectured to be non-Abelian 
is the v = 5/2 state^'^. In table I we list some of 
the proposed edge states for the J^ = § state and their 
quasiparticle properties. The edge states we consider are 
the Moore-Read state^^'^^ also known as the Pfaffian, 
its particle-hole conjugate the Anti-Pfaffian^^'^^ and the 
(331)-state89. See also Ref. 44. Of these the (331)-state 
is an Abelian theory. The proposed edge theories consist 
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Figure 8. The modulating function of four candidate states 
for the f = 5/2 state. The proposed states and corresponding 
quasiparticles are listed in table f. The parameters used for 
this plot are Uc = 5 • 10^ [m/s], •;;„ = 1.4 • 10^ [m/s], T = 1 



[mK], a = 2.4 [^m] and b = 2.1 [/.im]. 



F. u = 7/3 state 

The next state we look at is the i' = 7/3 plateau^'^. 
The trial states and the corresponding quasiparticles with 
lowest algebraic decay are listed in table II. These trial 
states are the Abelian Laughlin^^ state at j^ = 2 -f 1/3, 
the particle-hole conjugate of the Read-Rezayi^"^ state at 
k — 4, and two Bonderson-Slingerland states^*. The BS 
states are formed through a hierarchical construction of 
a non-Abelian candidate state, in this case the Pfaffian 
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Table II. Properties of quasiparticles for different edge models 
of the i" ~ I state. The e/3 Laughlin quasiparticle is present 
in the Laughlin 1/ = 2 + 1/3 state and all other non-Abelian 
states. Data obtained from Ref. 44. 
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Table III. Properties of quasiparticles for different edge mod- 
els of the V = ^ state. The 2e/5 Laughlin quasiparticle is 
present in all the listed states. Data obtained from Ref. 44. 



and Anti-Pfaffian state. Figure 9 shows the modulating 
function for the proposed states. 

In addition to plotting the tunnelling current for a 
number trial states, figures 8 and 9 show the effect of 
different values of g^ and gc on the R function. The gen- 
eral rule is that a larger value of gi corresponds to a larger 
damping on the contributing frequency. In particular, a 
larger sum of (?„ -I- gc corresponds to an R function which 
decays more rapidly for increasing V . 
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Figure 9. The modulating function of four candidate states 
for the u — 7/3 state. The proposed states and corresponding 
quasiparticles are listed in table II. The parameters used for 
this plot are Vc ^ 5 ■ 10^ [m/s], v„ = 1.4 • 10^ [m/s], T = 1 
[mK], a = 2.4 [^m] and b = 2.1 [/im]. 



G. 



12/5 state 



The last plateau we discuss is'* at z/ = 12/5. There 
are numerical studies^""^^ each of which suggest a dif- 
ferent quantum Hall trial state for the i/ = 12/5 
plateau. The edge states we discuss here are the particle- 
hole conjugate of the Read-Rezayi state^'^'^'^ at fc = 
3, a Haldane-Halperin edge^^'^'^'^^, and a Bonderson- 
Slingerland state^'*'^'^. The corresponding quasiparticles 
with lowest algebraic decay are listed in table III. The 
modulating functions for these states are plotted in 10. 



From the plots on the ly ~ 5/2, u ~ 7/3 and u = 12/5 
we find empirically that the parameters gi control the 
amplitudes of the different oscillations present in the R 
function. These are the oscillations discussed in Sec- 
tion Vine. We find that a larger gi causes a relatively 
smaller amplitude of the corresponding oscillation. This 
empirical rule is supported by the discussion on dephas- 
ing in Section VII D. Here it was found that for a typical 
length or temperature scale the R function is exponen- 
tially suppressed as a function of increasing temperature 
or increasing circumference of the interferometer. These 

scales are partially determined by [J2i ^\ ■ Here we 
find empirically that also the relative amplitude of each 
oscillation is inversely related to the corresponding alge- 
braic decay. 
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Figure 10. The modulating function of four candidate states 
for the v = 7/3 state. The proposed states and corresponding 
quasiparticles are listed in table III. The parameters used for 
this plot are v^ ^ 5 ■ 10^ [m/s], ■;;„ = 1.4 • lO'^ [m/s], T = 1 
[mK], a = 2.4 [^m] and b = 2.1 [fim]. 
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Figure 11. Idea of the setup of an interferometer with a side 
gate. By applying a voltage on the side gate the electrons 
are repelled thereby deforming the edge of the quantum Hall 
liquid. As a function of the side-gate voltage the effective area 
of the interferometer and the length of the lower edge grow or 
shrink. This changes both the AB phase and the R function. 



IX. THE AHARONOV-BOHM EFFECT AND 
THE INTERFERENCE CURRENT 

A. Weak tunnelling and the AB phase 

The AB-phase is determined by the magnetic field 
strength B, the area of the interferometer and the quasi- 
particle charge. It is given by the number of unit flux 
quanta for a quasiparticle with charge Q piercing through 
the interferometer 



pi*/*Q 



where 



$ = 2ttB X Area 



h 
Qe 



(99) 



Here <& is 27r times the total number of flux quanta 
through the interferometer and ^q is a unit flux quan- 
tum for a quasiparticle with charge Qe. This expression 
only applies in the weak tunnelling limit, where quasi- 
particles with the smallest algebraic decay are the most 
relevant operators in the language of the Renormaliza- 
tion Group. In this limit the interferometer is said to be 
in the Aharonov-Bohm regime and throughout this work 
we assume this always applies. 

In contrast, in the strong tunnelling limit the tun- 
nelling current effectively pinches off the area within the 
interferometer, thereby forming a quantum dot. This 
is called Coulomb blockade^*'^'"'. In this limit electrons 
tunnelling between the quantum dot and the fluid outside 
the interferometer form the most relevant operators. The 
AB phase is no longer determined by expression (99), see 
e.g. Ref. 96 for the case of the integer QHE. 



B. Manipulating the AB phase through a side gate 

The AB phase is manipulated by either varying the 
magnetic fleld strength or deforming the effective area of 
the interferometer. We are interested in the latter case. 
In practice"^ ^'■^'^ the area is changed through a side-gate 
voltage. This setup is depicted in figure 11 with the side- 
gate voltage given by Vg (not to be confused with the 
voltage bias between the two edges, ojq). By charging the 



side gate the Coulomb interaction repels electrons inside 
the interferometer, effectively deforming the area of the 
quantum Hall fluid. If we ignore the interference effects 
due to the R function or quasiparticle braiding, then the 
current shows the following oscillating behaviour due to 
the AB phase 



If 



Io+Ioscy<COsi<^AB{Vg)/'^Q+d) (100) 



<^AB{Vg) = 2tt 



B 



h/Qe 



X ATea{Vg) 



This oscillating signal arises in the weak tunnelling limit. 
One typically assumes the change in area is linear with 
respect to the side-gate voltage, meaning Area(V^) ex Vg. 
The Coulomb interaction and localization effects can al- 
ter this behaviour and cause small, non-linear fluctua- 
tions as a function of the side-gate voltage^^. This is 
called the Coulomb dominated regime (not to be confused 
with Coulomb blockade). In this regime the edge and 
the area inside the interferometer readjust to keep the 
dot neutral. Quasiparticles still tunnel along the point 
contacts and the interference current is still visible, but 
the corresponding AB phase does not follow expression 
(100). We assume the interferometer is not Coulomb 
dominated and the change in area is linear with respect 
to the side-gate voltage. 

Recent experiments^'^''^^'^^ observe on the order of < 5 
full oscillations when the side-gate voltage is varied. This 
applies to the v = 5/2 state, with a magnetic field 
strength of B ^ 5.5 [T]. The area of the interferome- 
ter is estimated in the range ofO.l - 0.4 [/im^], depending 
on the exact geometry of the device. For an interferom- 
etry area of 0.15 [/im^] and a quasiparticle with charge 
Qe — e/4 this corresponds to a total of roughly 50 unit 
fiux quanta. If we assume the interferometer is in the 
AB dominated regime, then a generous estimate of the 
change in area is about < 10% when five full oscillations 
are observed. 



C. The interference current: combining the AB 
phase and the H'^° function 

When the side-gate is used to change the area of the 
interferometer, then almost inevitably the length of the 
edge between the two point contacts changes as well. 
This change in length causes interference effects through 
the modulating function Hf^°'^. Including this in the ex- 
pression for the interference current gives 



Ip. 



^osr. ^ RG 



e'^AB{V,)/^Q+^5JJU,od^^^^.^y^■^ (^Q^) 



The function Hff^ 



implicitly depends on the side-gate 



voltage Vg through the length of the lower edge, aiVg). 
Whether the change in Hf^°'^ as a function of Vg is sig- 
nificant is determined by the change in the length of the 
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Figure 12. Plot of estimated AB oscillations as a function of 
the varying edge length. This variation is caused by the side- 
gate voltage and we assume a linear relation between the area 
of the interferometer, the side-gate voltage and the length of 
the edge. The parameters used are Q = e/4, gc ~ j^^ 9^ ~ \^ 
Uc = 8 • 10^ [m/s], «„ = 3 • 10^ [m/s], h = 2.5 [/xm], a = 2.25- 
2.75 [^J.m\, V^ = 50 [^V], T = [K]. The plot is "weak" in 
the sense that the modulating function does not change much 
over the plotted range. 



Figure 13. The same plot as in figure 12, but with different 
velocities, namely v^. = 9- lO'^ [m/s], Vn — &■ 10^ [m/s] and the 
range over which the edge length is varied is larger, a = 2.0- 
3.0 [/^m]. These slower velocities and larger range lead to a 
Hf^°'^ function which varies significantly more than that of 
figure 12. 



D. Frequency analysis of interference current 



edge, the velocity of the edge modes and the voltage bias 
between the two edges loq. 

For instance, in the experiment of Ref. 35 the quantum 
Hall fluid inside the interferometer is cigar-shaped with 
the ends of the cigar corresponding to the point contacts. 
We can picture the scenario in which the side-gate voltage 
deforms the lower edge uniformly, such that a 5% change 
in the area of the interferometer is accompanied with 
relatively negligible change in the length of the edge. In 
this scenario the function H"^°'^{luq; Vg) is approximately 
constant as a function of Vg. 

The other possibility is that the change in a{Vg) is not 
small. The device used in the experiment of Ref. 33 has 
a circular shape, and it is possible that the change in 
edge length is relatively larger than that of Ref. 35. It 
then depends on the remaining parameters, the velocity 
and voltage bias, if the change in Hf-°'^{ujQ; Vg) is large 
enough to be observable. 

In figure 12 and 13 we have plotted these two scenar- 
ios. Figure 12 is the "weak" case in which the function 
^mod j-gu^ajns largely constant while Vg is varied. In 
the lower panel of this figure the function Hf^°'^ causes 
a small modulation of the total interference signal. The 
interference due to a varying edge length is difficult to 
observe through measurement of this signal. Figure 13 
shows the "strong" case where the change in Hf^°'^ is 
much larger. These plots differ in the values used for the 
velocities and edge lengths, keeping all other parameters 
fixed. The frequencies of the oscillations at which Hf,-°'^ 



In plotting the figures 12 and 13 we assume a lin- 
ear relation between the area and the side gate voltage. 
Area oc Vg, and the length of the edge and side-gate volt- 
age, a{Vg) ex Vg. Under this assumption the interference 
due to the AB effect oscillates at some frequency with 
respect to the varying edge length a{Vg). We denote this 
frequency by cJ)ab, 



gi*/$Q _ g2iri0Af 



■a{V, 



(102) 



In other words, (pAB corresponds to the frequency of the 
oscillations appearing in the upper panels of figures 12 
and 13. Fixing the proportionality constant between the 
change in area and the change in edge length equal to 
Ci, i.e. AArea(Vg) = Ci x Aa(Fg), then (J)ab = Ci x 
The proportionality constant depends on the exact 



h/Qe 

details of the interferometric device, and the change of 
both area and edge length is performed through the side- 
gate voltage Vg. The charge of the quasiparticle in the 
fractional regime can then obtained by looking at the 
ratio of this frequency compared to that in the integer 
regime where Q — 1, 



bAB{v = fractional) 
(j^Asiv = integer) 



(103) 



varies are given by Eq. (98). 



In figure 6 we showed that the as a function of a vary- 
ing edge length the modulating function Hf^°'^ oscillates 

with frequencies "^r- ^ similar analysis shows that the 
combined signal of the AB oscillations and the modulat- 
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Figure 14. The four figures on the left plot the interference 
current Re[e"*'^^"ii"™°'*] as a function of varying the edge 
length a = [2.0 — 3.0] [/im] for an applied voltage bias of 
Vbias = 10, 20, 30 and 40 [^V]. The voltage bias Vbias should 
not be confused with the side-gate voltage. The figures on the 
right are the corresponding Fourier transforms. The peaks 
correspond to Eq. (104). The remaining parameters are Qc = 
g„ = 1/8, T = [K], Qe = e/, 6 = 2.5 [^m], «„ = 5 • lO" 
[m/s], ^c = 9 ■ 10" [m/s]. Finally, (J,ab = Ci x ^ =4.11 

[^m^^] with Ci = 1.0 • 10"^ [Aim] and B = 6.8 [T]. 



ing function oscillates at three frequencies, given by 



frequency peaks : 



'}ABi9AB 



QeV 

Vnh ' 



'JAB 



QeV 

Vch 



(104) 



These frequencies correspond to the signals appearing 
in the lower panels of figures 12 and 13. In particular 
the "pure" AB oscillations corresponding to ((iab remain 
present and the quasiparticle charge can be measured 
through formula (103) even if the oscillations in Hf^°'^ 

are strong. The shifted peaks (pAB + ^^ provide an 
experimental probe of the velocity of the edge modes. 
Numerical estimates^^ indicate that Vc > Wn, meaning 
the largest frequency corresponds to the velocity of the 
neutral mode. 

If the terms ^^ in (104) are small compared to <J)ab, 
then the frequency peaks overlap in a Fourier analysis 
and become indiscernible. To enhance the visibility of 
the different peaks we can either reduce 4>ab, increase the 
range over which the edge length is changed or increase 
^^^ . The frequency (J>ab and the variation in edge length 
are both determined by the geometric properties of the 
interferometric device and the side-gate voltage. 

Increasing "^^ can be accomplished by injecting a 
larger current into the system which is equivalent to in- 
creasing the voltage bias luq. The effect of a larger cur- 
rent on the frequency spectrum is demonstrated in figure 
(14), where the applied voltage bias is increased from 10 
[fiV] to 40 [fiY] in steps of 10 [fiV], while keeping all 



other parameters fixed. For a voltage bias of 10 [fiV] the 
frequency peaks merge and are indistinguishable. This 
is due to the relative magnitude of ^^ and ipAB- At 
a voltage bias of 40 [/iV] three frequency peaks emerge 
corresponding to the frequencies (104). 



X. DISCUSSION 

We have calculated the tunnelling current through a 
Fabry-Perot fractional quantum Hall interferometer in 
linear response theory for a broad class of edge theo- 
ries. Our main result is an expression for the tunnelling 
current in terms of Carlson's R function at finite tem- 
peratures and in terms of the confluent Lauricella hyper- 
geometric function at zero temperature. This expression 
arises as the tunnelling current is related to the Fourier 
transform of the quasiparticle propagators. In the confor- 
mal limit these propagators have a universal form, which 
is the reason behind the generality of our result. 

Our result applies to both Abelian and non-Abelian 
edge theories with an arbitrary number of edge modes 
- the neutral and charged degrees of freedom - each of 
which is characterized by its own edge velocity and chi- 
rality. In addition our result is applicable to interferom- 
eters with different edge lengths between the point con- 
tacts and our result can be straightforwardly extended 
to include more than two point contacts as explained in 
Section IV. 

We have implemented a numerical scheme to calculate 
Carlson's R function and the confluent Lauricella hyper- 
geometric function, and the corresponding interference 
current using a series representation. This numerical 
scheme is written in NumPy and Fortran and publicly 
available^^. We are also making available all the code 
that reproduce the plots in this work. 

The interference in the tunnelling current is attributed 
to the Aharonov-Bohm phase, the dynamical interfer- 
ence induced by the voltage bias between the edges and 
the statistical properties of the quasiparticles. Recent 
expcrimcnts^'^''^^ measure the Aharonov-Bohm phase by 
deforming the area inside the interferometer through a 
plunger gate. This setup also changes the edge length 
between the point contacts which induces interference ef- 
fects through the dynamical interference. We show that 
the total interference results in oscillations in the tun- 
nelling current as a function of the edge length. We have 
determined the frequency of these oscillations in terms 
of the edge velocities and the source-drain voltage, i.e. 
Eq. (104). These frequencies can be used to measure the 
edge velocities. 

The visibility of the frequency peaks depends among 
other things on the geometry of the interferometer and 
the range over which the length of the edge is varied. If 
there are many AB oscillations, while the change in edge 
length is small then the dynamical interference effects are 
hardly discernible from the AB oscillations. It is possi- 
ble that the change in edge length of current interfero- 
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metric devices is negligible and the interference effects 
we describe are indeed not measurable. In this case our 
proposed experiment requires an alteration of the inter- 
ferometer, using for instance a different geometry. 

Another way to increase the visibility of the frequency 
peaks is to increase the strength of the injected current. 
The frequencies are directly proportional to the source- 
drain voltage bias. The frequency of the AB oscillations 
is independent of the source-drain voltage, while the fre- 
quencies of the oscillations due to dynamical interference 
increase with larger source-drain voltages. 

Alternatively, this dependency can be used to check 
if the dynamical oscillations play a role in experiments 
which focus on the AB interference. If the effect of 
dynamical oscillations play a role in experiments which 
measure AB oscillations, then our results predict that 
this becomes apparent by running the experiment multi- 
ple times at different source-drain voltages. 
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Appendix A: Fourier transform of propagators 

We require the Fourier transform of the two-point 
propagator. We start with the finite-temperature case 
and no spatial dependence. 



Pg{t) = {-kTY sin ['kT{5 + it)\ 



(Al) 



Here ^ > is an infinitesimal integral regulator which is 
taken to zero in the end. To compute its Fourier trans- 
form Pg{Lo) = J dte'^^^Pgit) we follow Ref. 98 and substi- 
tute 5 + it = i-. 



'27tT 



^. This leads to 



Pg(t) = {irTycosh 



(f) 



= (27rT)^e"S^/^(l + e-^)" 



(A2) 



With this substitution the limits of the contour are 
±oo -|- i(7r — S). The contour is deformed so that it runs 
over the real line of x, which can be done provided there 
are no singularities that prevent this deformation. The 
function cosh(a;/2) is zero at the points a;„ = {2n+ l)7ri, 
for n integer. It is therefore the presence of the integral 
regulator (5, which allows for the deformation. 

After substitution the resulting integral is an integral 
representation of the Euler beta function^". We have 



PM = 



Q2T 



(27rr) 



i-s 



-(f-»2w)=^(l + e-^) ^ dx 



g2T 



(27rr) 



i-s 



B 



+ i 



2ttT 



2ttT 



(A3) 



Here we have taken (5 — >■ in the final result. 

Through a similar manipulation the Fourier transforms 
of products of two-point propagators with unequal ar- 
guments can be obtained. This results in the Fourier 
transform for G^, Eq. (65). We first note that with the 
substitution 5 + it — i 2^ -I- 2^ we have 

Pg {t + = (27rT)^ e-^^^^e-i-AgiO (A4) 

where Ag{0 = (l + e-^-^^e^")"^ . 

When applying this substitution to the Fourier transform 
of the product of four propagators (setting g = gn + gc) 
we obtain 

dt[e-*P,,(t-f^)P,,(t-A) 

e#(2^r)2^-ie^^('-")(S+''S:) 

/oo 
d:c[e-^9~^^)-AgAf-JAgA~^J 
-CSO 

X AgJ^j^JAgA-^l)] . (A5) 

The resulting integral is an integral definition of Carlson's 
R function^^, see Eq. (B3). This function is a scaling 
function and is closely related to the Lauricella hyperge- 
ometric function^®'^^ Ff^ . This Lauricella function is a 
multivariable generalization of the Gauss hypergeometric 
function of one variable and the Appell hypergeometric 
function of two variables^^. For our purposes it is conve- 
nient to use the R function to represent our main result, 
although the two representations are interchangeable, see 
Eq. (B9). 

Applying the integral representation (B3) gives for the 
integral 



(A5)=e-^(''-''Kft+')e)p2g(c.) 
xR(^g 



2t:T 



{gc,gc,gn,gn}; 



-2ttT- 



2-kT- 



-?)27rT- 



r]2TTT- 



(A6) 



where Pg{uj) is given by (A3). The resulting R function 
is Carlson's R function. Note that the order in which the 
parameters appear is important. 



Appendix B: Carlson's R function 

1. Main properties of the R function 

We first introduce a notation. We define Qn as the 
ordered set with n elements given by 



Gn = {gi,...,5n} 



(Bl) 
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and we set 






(B2) 



In the main text we usually work with the case where n — 
4 and the ordered set corresponds to Qn — {Qc, Qc, gmdn}, 
and 7 = '2{gh+gn) — 2g. Carlson's R function is treated 
in Rcf. 57 and is defined through the integral representa- 
tion 



R{a;gn;{zi}) = 



1 



i? (a,7 — a) 



dx . (B3) 



Here, B(x, y) is the Euler beta function, and {zi\ is the 
ordered set {zi}^^^ = {zi, . . . 2:„}. 

We require Re [a] > and Re [7 — a] > for conver- 
gence of the integral. Furthermore, we take the z^'s to 
be real and positive. The R function is symmetric under 
the simultaneous interchange of gt -f-^- gj and Zi ^ Zj . In 
the text the Zi correspond to the exponentials e ^ "^^ ^i . 

The R function is a scaling function, i.e. it is homoge- 
neous. This follows directly from the integral definition 
(B3) 



R{a;gn;{zi,...Zn}) = 

A"-R(a;^„;{Azi,...,Az„}) 

We also have the Euler-type transformation 
R{a;gn;zi,. ..z„) = 



(B4) 



n 

[n 



i?(7-Q;t;„;z^ \ 



(B5) 



For some special values the R function with n arguments 
reduces to one with m < n arguments. For instance 



Ria;g„;{zi,...Zk,z,. 
R{a;{gi,. 



,4) = 

■ ,gk,g};{zi,---Zk,z}) (B6) 



where g = ^fe+i + ■ 



5n • We also have the case 



B{a, 7 - a)R {a; C/„; {z^, . . . Zfe, 0, . . . , 0}) = 

B{a, 7 - a - g)R {a; Q^]{zi,... z^}) (B7) 

The R function is closely related to the Lauricella hyper- 
geometric function^^"^^. We define the Lauricella func- 
tion through its series representation 



^i)" (a;5n;7;{i -wi,...\ - w„}) 
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(B8) 



where {a)m = T[a + m]/T[a] is the Pochhammer symbol 
and we require |1 — w^j < 1 and arg(l — Wi) > for 
convergence of the series. 

To demonstrate the relation between the two functions 
we define z„ = max(zi, . . . , z„) as the largest parameter 
of the Zi's. Because of the identity (B3) we can always set 
this parameter to be the last argument of the R function. 
Furthermore, we will demand Zi ^ Zj for i ^ j, which can 
always be accomplished through the reduction property 
(B6). The relation between R and Fd is given by 



Rioi;gn;{zi,...Zn}) = 



1-^^^}) (B9) 



The arguments of the Lauricella function all satisfy 
|1 — Zi/zn\ < 1 and arg(l — Zi/zn) — meaning we have 
convergence of the series (B8). 

The Lauricella hypergeometric function is a general- 
ization of the single- variable Gauss hypergeometric func- 
tion, denoted by 2-^1, and the two- variable Appell hyper- 
geometric function, Fi. We have 

i?(a;{5i};{zi})-zr" (BIO) 

R{a;{gi,g2};{zi,z2}) = 

Z2-"2i^i(«;{gi};gi-fg2;{l-f}) 

R {a; {gi, g2, gs}; {zi, Z2, Z3}) = 

Z3-"i^i(a;{5i,52};ffi+ff2+53;{l-t,l-ft})- 
Here we assume Z3 > Z2 > zi. 



2. High temperature behaviour 

Consider again the expression for Gfj{uj), Eq. (G5). 
This is proportional to the integral 

/OO ^^ 

dte''^'l[Pg^{t + ^,) (BU) 

where the £,i correspond to the energy scales set by the 
velocity and edge lengths, £, ~ ^Vif: and Pg{t) is given 
by Eq. (Al). We are interested in the behavior for this 
function when T grows large. For this we substitute 
5 + it ^ ix + ^ , which gives 

Pg. (t + e.) = (TTTy cosh (ttT [x + C^)^'^ (B12) 
and the integral becomes 

/OO 
dxe"^"" Y[ cosh (ttT {x + ^,))"^- 

To be consistent with the main text we set ^^ gi = 2g. 
We split the integral into two domains, and pull out an 
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exponential from the cosh function. This gives 
e 2T 



r°° e 

Jo 

J —oo 



27rT 



= -^TC,\-s 



)-»' (B13) 



Consider the first integral. We perform an integration by 
parts, and obtain a boundary term and a remainder, 



Jo 

^ 2TTTg - iojl^^^^ 



-27rT(x+f )N-5i 



T«, 



e 



-TrT^i 



+ / da;e-(2^^»-^'^)^/(x) (B14) 
where 

X / 9i J— 

^—^ QTvTii _|_ g-27rT(x+^) 

We can estimate an upper bound for the remainder term. 
For this we note that f{x) is positive on the integration 
domain and bounded by 

f{x) < AnTg [|(e''^«' + g-'^^?' )~''' , 2: e [0, 00) 

i 

This gives an upper bound on the remainder given by 
1 



2TTTg — ioj 
AnTg 



dxe-^^'''^3-"^^'=f{x) 



< 



Y[{e^^^^ + e-^'^'^r 



, 2-Ti^)-9' / dree- (2-^9-'")^ 
27rTg - iw -"-.^ Jq 

AnTg 
{2TTTg - iuY 



J](e-Te. +e-T5.)-5. (B^5) 



The product also appears in the expression for the bound- 
ary term in Eq. (B14). This product therefore determines 
the asymptotic behavior of the boundary term in the high 
temperature limit, and also acts as an upper bound on 
the remainder term. A similar analysis can be applied 
to the second integral in Eq. (B13). It follows that the 
asymptotic behavior of the integral T in the high tem- 
perature limit is given by 



(27rr) 



2g-2^-7rTi:, IC.Iff. 



(B16) 



The factor (27rr)^^~^ is the high temperature behavior 
of expression (A3). This shows that the high tempera- 
ture behavior of the modulating function is given by the 
exponential exp(-7rrX]i ICi|ff»)- 



3. Computing the 7? function 

For n = 1 and n = 2 the R function reduces to the 
Gauss and Appell hypergeometric functions respectively 
for which various efficient numerical implementations ex- 
ist. For n > 3 no numerical implementation is available 
and we can either perform numerical integration or com- 
pute the expansion (B8) to some finite order. Numerical 
integration of the integral (B3) takes into account the 
Beta function as well, which is why we use the series 
expansion instead. We will follow Ref. 84 to cast this 
series expansion into a more tractable form suitable for 
a numerical implementation. 

The main result of Ref. 84 is that the multivariate Tay- 
lor expansion (B8) can be written as the single summa- 
tion 

F}^' {a;Qn--,l\{l -wi,...l- Wn}) 
00 , ^ 

Here (a)„ — T[a + n]/r[a] is the Pochhammer symbol 
and Am is the cycle index (of the symmetric group Sm) of 
the variables tj. Defining the variables tj (j = 1, . . . , m) 



^3 =^9^{^ 



Wi 



(BIB) 



4=1 



then the cycle index A^ of this set {ij}"Li is given by 



^^7n i.^1 5 • ■ • ; ^m) 



E 



k\ ,...,krri 

ki-\-2k2 hmkjj-^—m 



-^J- k,i \j 



(B19) 



The summation over the fc^'s (B19) is constrained by 
^- j^ jfcj = 7n, which makes its computation for large 
TO rather involved. It's more efficient to use an iterative 
approach, as A^ can be expressed in terms of {A„}„<m. 
Defining Aq = 1 we have for m>\ 

Am(^l, • ■ • ,^m) = — 2^^J^™-j(*l' • • • '^rri-j) (B20) 

Let us also give the corresponding expansion for the R- 
function. For that we again assume z„ is the largest 
argument of the function. Then 



R {a; {gi, . . . , gn}; {zi, . . . , Zn}) = 



(B21) 



where 7 — J27=i 9i ^^'^ ^n — max(zi, . . . , z„) as before. 
The Tj {i — I, . . . ,n ~ I) are given by 



n-l 



Tj = ^giil- z,/z„y 



(B22) 



i=l 
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This algorithm is due to Laarhoven and Kalker®^. 

In the main text the R function which enters 
the expression for the interference term is a multi- 
variate expansion in terms of the scales 1 — -£i- = 
1 — exp(— 27rT(— — — )) < 1. For large temperature 

scales (> 15 mK) the arguments approach the radius of 
convergence, (1 — Zi/z„) < 1, and the rate of convergence 
of the series becomes extremely slow, especially when the 
frequency ljq becomes large as well. This requires a very 
large number of terms in the expansion, which becomes 
problematic since the algorithm for A^ scales as order 
0{N^) with N the number of terms in the series. In this 
regime numerical integration does not seem to be an al- 
ternative, as the standard integration schemes suffer from 
slow convergence as well. 

The situation is somewhat improved by using a series 
acceleration. We have chosen a series acceleration via 
the Combined Nonlinear- Condensation Transformation 
(CNCT) as outhned in Ref. 88. The algorithm works in 
two steps. First, the (largely monotone) series (B17) is 
transformed into an alternating series via a Van Wijn- 
gaarden transformation. Alternating series are known to 
converge better using a series acceleration. Second, this 
alternating series is accelerated via a nonlinear sequence 
transformation. For our purposes we have chosen Levin's 
u transformation*^, although other choices yield similar 
results. 

The advantage of the CNCT method is that only a 
handful of terms of the original series are needed to ob- 
tain a high precision estimate of the series. This method 
significantly improves the rate of convergence of many 
series**. However, the method requires the capability to 
compute "random" terms in the series (B21). To be spe- 
cific, to perform the Van Wijngaarden transformation we 
require the terms t^MAm with M — 2''{j + 1) — 1 and j 
and k integers, see Ref. 88. Typically we need all terms 
with j, fc < 30 for a decent precision in the final answer. 
But note that the index M grows exponentially. This 
is problematic, because our algorithm is designed to de- 
termine Am iteratively and this iteration process grows 
as 0{N'^). The CNCT method and similar acceleration 
methods therefore do not fully resolve the issue of slow 
convergence. To avoid this problem our plots are per- 
formed at low temperature {T — 1 [mK] or T = [niK] ) . 

A second problem that arises is a lack of precision in 
the terms computed. We found that the typical dou- 
ble floating point accuracy can lead to problems when 
evaluating the series for large ujq (> 100 [mK]) and val- 
ues of the velocities and distance scales as mentioned 
in Section VIII. This issue is resolved by making use 
of high-precision floating point accuracy^^. The down- 
side to this is that the computation of a large number of 
terms is extremely slow. In particular, we cannot simul- 
taneously make use of the CNCT algorithm and high- 
precision floating point accuracy. 

We have implemented this algorithm through a combi- 
nation of NumPy^"" and Fortran, making use of F2PY^"^. 



In some cases we also made use of the high-precision 
floating-point arithmetic package mpmath^^. All plots are 
generated using matplotlib^"^. 



Appendix C: Zero temperature case 

In the zero temperature case the KMS relation of the 
G correlators, see Eq. (39), no longer applies. Within our 
approximation we do have the relation GfJt) = Gf^{—t). 
The expression of the tunnelling current at zero temper- 
ature is therefore given by 

N 

Ib{uJq) = ge(^ |r,|2 [G>{uq) - G<{u:q)\ + 



N 



2^|r,r*|Re[e 



i^ij+iccij l/^> 



[G>{ujQ)~G<iu;Q)]] 



i<j 



(CI) 



The analysis of the G^ correlator is the same as in the 
flnite temperature case, with the exception that we use 
the zero temperature expression of the propagator Pg (t) . 
In particular (65) still applies, but with the propagator 
given by 



m) = 



1 



{S + it)9 ' 

The expression for G^ and G^ then boils down to 



(C2) 



Gf-ji^) - Gfiiuj) = a^^,cAij{x)v, 



-2g,^-2g„ 



die" 



„/ o. , „ , &,„. a^„. b . 

PgAt + -)PgAt - -)Pa„ {t + V- Pg„ {t - V- h 

Vc Vc Vn Vn 

„/ a ^ ^ , b . ^ , a . ^ , b . 

PgA~i+-)PgA-t )P9A-i + V-)PgA~t-V-) 

Vc Vc Vn Vn 

(C3) 

We have not found a reference or method to treat this 
Fourier transform directly. It can be treated for the spe- 
cial case of a symmetric interferometer and a single edge 
mode, where Vc — «« and a = b. This special case is 
treated in Appendix C 1. Alternatively, we can start with 
the expression for the tunnelling current of the flnite tem- 
perature case and take the zero temperature limit. This 
approach allows for more general values of the physical 
parameters and is performed in Section C 2. Finally, we 
suggest in Section C3 a solution to the integral (C3), 
obtained by taking the zero temperature limit from the 
flnite temperature expression. 

As in the finite temperature case we find for the tun- 
nelling current 



-^b(wq) = 



Qe 



2Sc 2g, 
Vc Vn 



-avac|r(a;Q)|2/2g(wQ)sgn(wQ) (C4) 
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with l2g given by (C7) and the effective tunnelhng am- 
plitude equals 



Af 



i=l 

N 

2^ |r,r*|Re [A,(x)e'*'^+^"-i?r'(^Q)] (^5) 



i<j 



The modulating function Hf^ {uq) is given by (C13) in 
the symmetric interferometer case with a single mode, 
and by (C19) in the more general case. 



1. Zero temperature: tunnelling current for a 
single mode by computation of the Fourier transform 

We start with the Fourier transform of the correla- 
tor Gii{uj). This corresponds to the tunnelling current 
through a single point contact, see Eq. (CI). We require 
the Fourier transform of the propagator Pg{t), which is 
given by^'^ 



where 






Ig{L0)Q{L0 



w 



9-1 



(C6) 

(C7) 



and 8(a;) is the step function. Then 

G>{uj) - G<{uj) = a.aci'-'^^t'^'^"/2g(w)sgn(w) . (C8) 

For the expression of the interference term we set 
V — Vc ~ Vn and a = b. The required integral is (see 
Eq. (C3)) 



Kgiu;) 



V V 



a.. 
-t+-)] 

V 



(C9) 



We consider the separate cases where g <1 and g > \- 
The two cases overlap, and we find a single expression 
applicable for all values of 5. For <? < 1 the integral regu- 
lator is not required, so we set 5 — i). With some careful 
manipulations of the fractional powers of i we obtain 

sin(|a;|t) 



K„<,d^) = 4 sin(7r.g)sgn(w) 



dt 



{t'-ilY) 



9L\2\g 



1 /M^^^-9 

L-^ 2^ V 2w 



X J,_i(4^)/2,(^)sgn(^). (CIO) 



The function Jg{x) is the Bessel function of the first kind. 
The integral is found in Ref. 60. For the case of g > 5 we 
need an integral representation of the confluent hyperge- 
ometric function i.F'i 



e-'"iFi(5; 2g- (7 - P)uj)hg{i^)Q{i^) . (Cll) 



This applies when Re [/3], Re [7] > and Re[g\ > \. With 
this integral representation we find 



i^,>.(c.) = e-^He^^ 



XiF] 



(5; 25; 



-2V- 



l2giLo)sgn{Lu) . (C12) 



For these specific parameters the confluent hypergeomet- 
ric function iFi reduces to the Bessel function of the first 
kind^^o 



li^i 



(3; 25; 



-2i- 



\uj\a 



9 



f\uj\a\h-9 f\uj\a\ 



Therefore both cases (5 > 5 and g < 1) match in the limit 
of (5 — > and expression (CIO) extends to all values of g > 
0. Finally, we have for the zero temperature expression 
of the modulating function of a symmetric interferometer 

HrM = rb + l](!^)*-Vi(it)- (C13) 



2. Obtaining the zero temperature limit from the 
finite temperature expression 

The more general case in which we consider multiple 
modes with different edge velocities involves a more com- 
plicated Fourier transform which we are not able to de- 
termine directly. Instead, we use the result for finite tem- 
peratures and take the limit of T | 0. 

We require the zero temperature limit of the modulat- 
ing function, Hf^°'^ see Eq (69). To perform this limit we 
make use of the series representation of the R function, 
Eq. (B8). This gives 






. ^27rTxi 



, e 
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00 



lime-(2-^"-'-)-" y... > 

mi— niTi-i— 
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YT^Zi m. 



n-1 

n 
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(CM) 



The Xi's correspond to the (real valued) energy scales 
associated with the edge modes, i.e. ±— and so on. We 
assume a;„ > Xi for all i and we write x^ — Xi = Xn.i > 0. 
The limit is determined term-by-term. We first note the 
approximation 



1 



'2TrTx„ 



■')"" = (2^T)'"'x:::; + 



(C15) 



The dots are of higher order in T. Combining this with 



the (a - 



-27rT 



)m term, where m = X]i=i "^«i ^^ obtain 
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for the zero temperature limit 



follows from the series representation. And so 
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m—l 
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(C16) 



The higher order corrections of (C15) vanish in this limit. 
Plugging this back into (C14) gives 



lim i? (a - i 2^ ; 5„ ; {e'^^^^ , . . . , e^^^^" }) 
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with in — X]r=i "^j- ^^^ resulting series is called the 
confluent Lauricella hypergeometric function^^ 



lim R { q — i — — 

Tio V 2ttT 



;{5,5};{e^-^%e-^n 
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(C20) 



This matches with the result (CIO). 

The series expansion of the confluent Lauricella func- 
tion (C18) is of the same form as the non-confluent Lau- 
ricella function, (B8). The same combinatoric trick as 
explained in Appendix B 3 can be used to rewrite this 
multivariable series as a single expansion in terms of cy- 
cle indices, see Section B3. This expansion is given by 

4'"\Gn;r,{wu...,Wn})= (C21) 

We flnd that the convergence of the confluent series is 
much better than the non-confluent (finite temperature) 
case. In general, we do not require as many terms in 
the series. However, for the physical values of the veloc- 
ity, distance and voltage used in the main text we flnd 
that double floating precision is still not sufficient and we 
require high-precision fioating point numbers^^. 



'^2'\Sn-ii;{wi,...,Wn}) 



E 



00 

E 



7^0 m„=0^'^^^i' 



ft i^wT' 



.1=1 



(C18) 



This series is a multivariable generalization of the con- 
fiuent hypergeometric $2 function"'. The expression 
for the modulating function is 



^2 iSn-i;i;{~ii^Xn,i, 



-iUJX„,n-l}) . (C19) 



Here we recall that the Xi correspond to all combinations 
oirji— and —rji — , the parameter Xn satisfies x„ > Xi for 
i < n and Xn^i = x„— x^ > 0. As a sanity check we look at 
the case treated in Appendix C 1, which corresponds to 
the symmetric interferometer and a single channel. The 
confiuent Lauricella function reduces to the confluent hy- 
pergeometric function, $2 (^: c; a;) — iFi{b,c;x), which 



3. Zero temperature: multiple modes 

We have obtained the general expression for the zero 
temperature case by taking the zero temperature limit 
of the finite temperature expression. The same result 
can also be obtained by taking the Fourier transform 
of the zero-temperature expression for the G^ correla- 
tors. Since these calculations must produce the same 
answer we obtain the following integral representation of 
the confluent Lauricella hypergeometric function. With 
Pg{t) = {S + it)~s we have 

/oo 
dte^-^[Y[Pg^it + x,)-l[P,,{-t + Xk)] = 
3 k 

/2g(w)sgn(w)x 

e'"""$^"''^(a„-i;7;-*wa;„,i,--- ,-iujx„,n-i) ■ (C22) 



Here a;,, 



Xi > for all i < n, all gi > and 6 



is taken to zero in the end. The function Ig is given by 

(C7) and $2 ^ i^ ^^'^ confluent Lauricella hypergeomet- 
ric function, which has the series representation (C18). 
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